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Abstract: We show that particle-hole pairing is realized in the background of a charged 
black hole in magnetic field. The pairing instability occurs for sufficiently large fermion 
charges, which correspond to the Fermi liquid regime. The critical temperature for Fermi 
liquids is proportional to the magnetic field and vanishes as we approach the non-Fermi 
liquid state. The pairing order parameter leads to a relative shift of the Fermi surfaces 
corresponding to the bulk fermions with spin up and down. The value of the shift in Fermi 
momentum kp and the critical temperature Tc are proportional to the effective density of 
states at the Fermi surface. Our one-loop calculations provide a dual description of the 
magnetic catalysis for the lowest Landau level in graphene. This analyses may be relevant 
for the antiferromagnetic behavior in the cuprate superconductors and for the chiral spirals 
in the chiral magnetic effect. 

We also discuss thermodynamic and transport properties of a system at the boundary at 
zero magnetic field. The scaling behavior of the specific heat is c ~ T for Fermi liquid 
and c ~ T^'^ for non-Fermi liquid, while the behavior of the DC conductivity is the same 
T-^" in both cases. While it can be difficult to extract transport and hydrodynamic 
from the lattice, the AdS/CFT approach provides a robust frame for nonperturbative 
calculation of these properties. 
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1. Introduction 

Particle- hole pairing appears in different contexts in condensed matter physics. We con- 
sider here magnetic catalysis, i.e., generation of the T-odd mass parameter in the presence 
of a magnetic field. It is a well established phenomenon in (2 -|- l)-dimensions and it is 
believed to explain the anomalous quantum Hall effect in graphene, i.e., the appearance 
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of the additional plateaus in the Hall conductivity axy for the lowest Landau level [1] (see 
[2] for the magnetic catalysis in (3 + l)-dimensions). Electron-hole pairing is responsi- 
ble for the spin density order parameter and the antiferromagnetic nature of the cuprate 
superconductors at half filling. Spin-density wave in the form of spin-orbit ordering can 
trigger the superconducting pairing, while both superconducting electron-electron and spin 
density wave electron-hole orders are essential to describe physics of Mott insulating and 
pseudogap phases. Recently, there was an interest to the particle-hole pairing in the form 
of chiral spirals in the context of the chiral magnetic effect and the quarkyonic matter [3] . 

These phenomena involve strongly coupled physics. We therefore use the AdS/CFT 
correspondence which is a powerful tool in understanding strongly coupled quantum field 
theories. It is formulated as a duality between classical gravitational theory in the anti-de 
Sitter [AdS) space and a strongly coupled conformal field theory (CFT) in the limit of large 
N and large 't Hooft coupling A defined on the boundary of the AdS space. Recently the 
AdS/CFT correspondence was applied to different phenomena which arise in the context 
of condensed matter systems [4, 5, 6, 7, 8, 9]. Many of the above studies were initiated 
by the original papers on a holographic superconductor [10, 11], the non- Fermi liquid 
behavior [12], and quantum phase transitions [13]. In particular, there have been significant 
developments in understanding the superconducting instability near a charged black hole. 
It was shown that charged black holes are unstable to forming hair, which means that 
a (free) charged (or neutral) scalar field develops a vacuum expectation value and breaks 
spontaneously the corresponding symmetry when put in the charged black hole background 
with asymptotic AdS geometry [10, 11]. This was linked to the Breitenlochner-Freedman 
instability that provides a gravitational mechanism for superconductivity: if the charge of 
the boson is sufficiently large compared to its mass it will condense. This mechanism does 
not give microscopic details behind the superconductivity like the BCS pairing does. It 
provides the evidence for the bosonic condensate and suggests a holographic mechanism 
for the superconductivity. Using the Cooper pair picture the critical temperature has been 
calculated in Ref.[14]. 

In describing the particle-hole pairing we follow the same route as used to address 
superconductivity [15]. Both graphene and the cuprate superconductors are systems at 
finite charge density in (2 -|- l)-dimensions. Therefore the gravity dual description is given 
by a charged black hole in (3 -|- l)-dimensional anti-de Sitter space-time AdS^. Strong 
coupling and large N limit of the boundary theory translates into a gravity theory at small 
curvature and low energy, which reduces to a universal sector of classical Einstein gravity 
plus matter fields. The global U{1) symmetry of the conformal field theory (CFT) with 
current is mapped to a U(l) local gauge symmetry with a gauge field Am in the AdS^. 
In the AdS^, Am is an actual (not background) ?7(1) field, which is dynamic. 

In this paper we add a four-Fermi contact interaction between the charged fermions. 
We choose the channel favoring the magnetic catalysis and look for the particle-hole in- 
stability that shows up when the one-loop effective action has negative modes. There is 
important difference between showing the instability for the bosonic field < ^ >^ and 
for fermions < tpip >^ 0. Calculation for bosons is classical in the black hole background, 
whereas for fermions it involves one loop computation. We use variational approach where 



- 2 - 



we utilize the formula for the one loop fermion determinant expressed through a sum over 
quasinormal modes of the black hole [16, 17]. As discussed in [16, 17], the quasinormal 
modes are given by the poles of the retarded fermion Green function. The structure of 
the poles for the retarded Green function has been obtained in [18] for various relative 
relations between the charge and the mass of the fermion. We also do a one-loop calcu- 
lation in the bulk to obtain a non-local in the radial direction Ginsburg-Landau action. 
The latter calculation involves bulk fermion propagators and a radial profile for the pairing 
order parameter. The idea of calculation follows the Ginsburg-Landau approach. 

In this paper we consider application of particle-hole pairing to the magnetic catalysis. 
Magnetic catalysis has been shown in (2 + 1) [1] and (3 + 1) [2] dimensional relativistic 
models. The general result is that a constant magnetic field leads to the generation of a 
fermion dynamical mass even at the weakest attractive interaction between fermions. The 
essence of the effect is that in the magnetic field the dimension of the system effectively 
reduces d — )• d — 2, i.e., to (0 + 1) and (1 + 1) dimensional systems, that favors the 
dynamics of the particle-hole pairing (therefore the name of magnetic catalysis) [1]. We 
can choose different forms of the four-Fermi interaction, that will generate a mass term 
for the fermions with needed symmetry properties. For simplicity, we choose a contact 
interaction Gj„,f (ViiF^r^-;/;)^ with the strength Gint written through the mass scale of the 
interaction Gint = l/-^i>it- We show that this interaction triggers the generation of the mass 
term A-^zF^r^^/^ which is odd both under time-reversal and parity transformations (see the 
representation of F matrices). Contrary to the Dirac mass term mipip^ the generation of the 
T-odd mass H^'ipiV'^T^ ijj does not break any symmetry, e.g. U{1) l x U{l)ji in the NJL model 
or spin (flavor) symmetry SU{2) (C/(2)+ x [/(2)_) in the case of graphene [1]. Because no 
symmetry is spontaneously broken by the T-odd mass, no gap opens in the spectrum. In 
this sense, the condensate A^zF^F^-i/^ is similar to the spin density wave < tp^aip >, where 
the former couples to the mass and the latter one to the chemical potential. 

The paper is organized in the following way. In section 2 we introduce the black hole 
geometry and consider the near horizon limit which is dual to the IR GFT. In section 
3 we perform the variational calculation for the particle- hole pairing order parameter. In 
section 4 we perform Ginsburg-Landau calculations in the bulk geometry, and cal- 

culate the critical temperature. In section 5 we consider thermodynamic and transport 
properties of a system on the boundary at zero magnetic field. Appendices contain solu- 
tion of the Dirac equation in magnetic field and calculation of the Landau levels in the 
AdS4 holography (Appendix A.l), calculation of the conformal dimension in the IR CFT3 
(Appendix A. 2), solution of the Dirac equation in the AdS2 and obtaining the IR CFTi 
conformal dimension (Appendix B.l), derivation of the AdS2 Green function (Appendix 
B.2), one-loop calculation in (2 + l)-dimensional field theory (Appendix C), calculation of 
the critical temperature in the AdS4^ (Appendix D). 

2. Dyonic black hole and infrared CFT 

We consider 3-dimensional conformal field theory (CFT) with global U{1) symmetry that 
has a gravity dual. At finite charge density and in the presence of magnetic field, the 
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system can be described by a dyonic black hole in 4-dimensional anti-de Sitter space-time, 
AdSi, with the current in the CFT mapped to a U{1) gauge field Am in AdS. 
The action for a vector field Am coupled to AdS 4, gravity can be written as 

where gp is an effective dimensionless gauge coupling and R is the curvature radius of 
AdS4. The equations of motion following from eq.(2.1) are solved by the geometry of the 
dyonic black hole, i.e., with both electric and magnetic charges. 



2 



dr^ 



ds" = gMNdx^dx^ = ^{-fdt^ + dx^) + (2.2) 

where the redshift factor, /, and the vector field Am reflect the fact that the system is at 
finite charge density and in the magnetic field, 

= (1 - ^) , A^ = -ny, (2.3) 

where we chose the Landau gauge; the chemical potential /.i and the magnetic field 7i are 
given by 

9fQ ^ gpH 
^=R^o'^ = ^- ^'-'^ 
Here rg is the horizon radius determined by the largest positive root of the redshift factor, 
f{ro) = 0, 

M = rl + (2.5) 

and the CFT is defined at the boundary r — t- 00. The geometry eqs.(2.2),(2.3) describes 
the boundary theory at a finite density, i.e., a system in the medium at chemical potential 
/i, with the charge, energy, and entropy densities given, respectively, by 

_o Q M _27rr|_ 

^~ K^R^gp' '~ K^R^' '~ R^- ^ ^^> 

The temperature of the system is identified with the Hawking temperature of the black 
hole, Th ~ |/'(ro)|/47r, 

3ro / Q^ + H^' 



Since Q and H have dimensions of [L]^, it is convenient to parametrize them as 

Q' = 3rt, Q' + H' = 3rt. (2.8) 
In terms of tq, r^, and r*,,, the expressions are 

f = 1 + ^ - ^0 +'^^**/^o 
At = n(l-'^), A:, = -n, (2.9) 
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with 

^ = V35F-^, ?^ = ^/35F^^5£f3. (2.10) 

The expressions for the charge, energy and entropy densities, and for the temperature are 
simphfied as 

r=AZ:^fl_!4V (2.11) 

In the first part of the paper we consider the zero temperature hmit, i.e., extremal 
black hole, 

r = ^ ro = r„, (2.12) 

which in original variables is + = Sfq. In the zero temperature limit, eq.(2.12), the 
redshift factor /, eq.(2.9), develops a double zero at the horizon 

f = Qt^I^ + ..., (2.13) 

As a result, near the horizon the AdS^ metric reduces to AdS2 x with the curvature 
radius of AdS2 given by 

R2 = ^R. (2.14) 

This is a very important property of the metric, which simplifies calculations. This metric 
reduction can be seen explicitly by considering the scaling limit 

_ Ri _T_ 

T r^^ — A^, t — 

A^OwithC, T finite, (2.15) 
then the metric eq.(2.2) describes a black hole in ^^5*2 x R? 

ds'' = ^{-dr^ + dC) + (2.16) 



^- = ^4?' ^x. = -^?/. (2.17) 



with 



Physically, the scaling limit eq.(2.15) with finite r corresponds to the long time limit of 
the original time coordinate t, which translates to the low frequency limit of the boundary 
theory 

^ 0, (2.18) 

A* 

where w is the frequency conjugate to t. (One can think of A as being a frequency lo). 
Near the AdS^ horizon, we expect that gravity of the AdS2 region of an extremal dyonic 
black hole is described by a CFTi dual. We refer to [18] for an account of the AdS2/CFTi 
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duality. In what follows we use the horizon of AdS2 region at ^ — ?■ oo (coefficient in front 
of dr vanishes at the horizon) and the infrared CFT [IR CFT) defined at the AdS2 
boundary, C = 0- The scaling picture eqs.(2.15),(2.16) suggests that in the low frequency 
limit, the 2-dimensional boundary theory is described by this IR CFT (which is a CFTi). 
The Green function for operator O in the boundary theory is obtained as a small frequency 
expansion and by a matching procedure of different regions along radial direction, and is 
expressed through the Green function of the IR CFT [18]. 

3. Variational calculations of the paring gap 

In this section we perform variational calculations of the particle-hole pairing gap in the 
bulk. The logic of calculations is the same as in a field theory, except for arising radial 
dependence of bulk quantities, e.g. for the gap parameter A(r) as opposed to the BCS with 
a constant gap. The radial profile is important to keep, since it insures convergence of radial 
integrals and for different bulk behavior characterizes different systems on the boundary. 
Our variational calculations in the bulk are possible due to the one-loop formula for an 
effective action obtained in [16, 17] and expressions for the poles of the fermion Green 
function obtained in [18]. 

3.1 Effective action for interacting fermions in a magnetic field 

We consider a spinor field "0 ™- the ^^5*4 of charge q and mass m, which is dual to an 
operator O in the boundary CFT^ of charge q and dimension 

A^ = ^+mR, (3.1) 

with mR > 2 and corresponding to the "stable" CFT. In the black hole geometry, eq.(2.2), 
the quadratic action for ■0 is written as 

So = i j d^x^/^ {i)T^^VM^ - mi^^f) , (3.2) 

where ■0 = ijj'^T-, and 

Vm = Om + ^u^abAiT"'' - iqAM, (3.3) 

with ujabM the spin connection, and F"* = ^[r'^, T'']; here M and a, b denote the bulk space- 
time and tangent space indices respectively, and /u, denote indices along the boundary 
directions, i.e. M = (r, /i). 

As discussed in the introduction, we can add to Sq the contact interacting part 

Sint = - j d^xV=5Gi„i(viir2r§0)(v;ir2r§^), (3.4) 

where Gint = ^/^inf ^int is a mass scale of the interaction. The representation for F 
matrices is given by eq.(A.15), and hat indices on F matrices always refer to tangent space 
indices. In this representation of F matrices, 

zF^F^ = 
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We also have iF^F^ = — F^F^F^. The form of interactfon eq.(3.4) is motivated by the form 
of the projectors eq.(A.37) which decouple ip into two components. 
We add the magnetic Zeeman splitting of the spin degeneracy 



Sb = J (Tx,/^ qni^ra-'i;, (3.6) 

where acts on spin indices. The resulting action is the following sum S = So + Sint + Sb- 
We solve the four- Fermi interaction, Sint, in the mean- field approximation by per- 
forming standard Hubbard-Stratonovich transformation. We introduce a composite order 
parameter 

A = 2Gint < i^iV^T^^ > (3.7) 
and decouple the interaction into a quadratic form 



Sint = J d'x^g - (A^/^iF^F^V' + /i-c.)J . (3.8) 

The order parameter A is T-odd, i.e. T{T^iT'^T^)T'^ = -(F^zF^F^), with 

and charge conjugation is fixed by CF* = F*". The exact form of the four- Fermi interaction 
is not important. With any interaction Gmj(V'F'i/')('0F'0) respecting the symmetries, where 
F stands for a collective combination of Gamma matrices, the gap given by eq.(D.4) can 
be generated. Its effective action is given by Seff = -iTi (inG-i + \{G^^G - 1)) which 
satisfies the stationarity condition (or gap equation) 5Seff/6G = [1]. 

To get an effective action for A, we intet grate out the fermion fields with the result 

S,ff = I d^x^g - ^TrlnG-i) , (3.10) 

where the full fermion propagator G{x,x') =< il){x)il}{x') > and its inverse is given by 

G-^{x,x') = T^^Vm -m- AiF^F^ ± qHVK (3.11) 

In the one-loop effective action eq.(3.13), the trace and the logarithm are taken in the 
functional sense. The coordinate x = {r, t, x} includes the radial r and (2 + 1) boundary 
directions {t, x} in AdS^. If we assume translational invariance along boundary directions, 
then the Fourier transform is given by 

G{x,x') = tY.J ^G(i.;„,A.,r,r')e— (3.12) 

where fermionic Matsubara frequency is Un = (2n + 1)ttT. We furthermore assume that 
the condensate is a function of only the radial coordinate in the AdS^, A{r,uj,k) = A(r). 
The effective action is given by 

.%„ = ^ - I E/ |^*'TVl„G-(i.„.*,.,.')] . (3.13) 
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In general, it is a difficult task to calculate one loop effective action in the bulk. Here we 
will use the method suggested in [16, 17] to calculate the free energy (determinants) in a 
black hole background as a sum over the quasinormal modes of the black hole. We will 
also use the recent analytical results [18] for the fermion quasinormal modes. 

3.2 Variational calculations of magnetic catalysis in a charged black hole ge- 
ometry AdS4^ 

To calculate one-loop fermion action eq.(3.13) we need to know eigenvalues of the Dirac 
equation eq.(A.42), which can be written symbolically for each mode as [16, 17] 

M{z,l)^ = X{z,l)^ (3.14) 

with z = iujn- The zero modes A(z*(/), /) = define as solutions the quasinormal frequencies 
z*(/). As was shown in [16, 17], the quasinormal frequencies of a wave equation in a 
black hole space-time are actually poles in the corresponding retarded Green function in 
the black hole background, where 

M{iiOn, l)G{iujn, I, r, r') = r^6{r, r'). (3.15) 

Indeed representing the Green function as a sum over eigenfunctions we have 

1 r+ 

Tr— -= drG{iujn,l,r,r), (3.16) 



' M{iu}n,l) 

which is the usual representation of G. This equation is shown to be true for general 
complex z = ioon and G satisfyes the ingoing boundary condition at w 7^ and regularity at 
a; = at the horizon [16, 17]. Therefore, as was shown in [16, 17], the fermion determinant 
given by a sum over the quasinormal frequencies of the black hole, i.e. when M{z,l) has 
a zero eigenvalue, is equivalent to the sum over the poles of the analytically continued 
to complex frequencies fermion Green function with ingoing boundary conditions at the 
horizon. This method is also used in color superconductivity [19]. 

Analytic results have been obtained for the Green function in the [18]. (Numeri- 

cally it has been obtained in [20].) A general form for the retarded Green function is given 
by [18, 21, 22] 

B+ + B^G'^ioj) 

^«("'^) = i7t:4:g^W)' ^ ^ 

where the ratio of numerator to denominator comes from expansion of the solution of 
Dirac equation at the boundary r — )■ 00, therefore the ratio B/A is the AdS^ retarded 
Green function, Gr. The coefficients A±,B± are expansions (rows) in small frequency 00. 
The low- frequency limit is equivalent to the expansion near the horizon at small oj, where 
metric reduces to AdS^ — )• AdS2 x R^- Therefore A^ {B+) and A^ (B^) components 
arise from expansion (and matching procedure) near the AdS2 boundary, that relates ratio 
A^/A- (B^/B-) via the IR Green function G^^ obtained in the AdS2 x R? calculation. 
Non-analytical frequency behavior of G is controlled by the IR GFT, G^^, while coeffi- 
cients A± and B± carry the UV information. In section 4, we find the coefficients A±, B± 
to the leading order in frequency, that requires solution of the Dirac equation in the AdS^. 
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It was found in [20] and [18], that the fermion Green function develops a sharp pick 
indicating the existence of the Fermi surface and quasiparticle poles. Expansion of the 
Green function near the Fermi surface is given by [18] 

GR{u,k) - 



S(a;, kp) = hvpG^^iuj, kp) = hvFc{kF)u'^''^F , (3.18) 

where we keep notations for constants introduced in [18]. Here k± = k — kp, hi, h, vp are 
governed by the UV physics, and were obtained numerically in [18] and here in section 4. 
Eq.(4.17) gives the rough structure for the boundary Green function. Further, in section 4, 
we use a more detailed description for Gr. As T — t- 0, Yi{ijj,kp) — ?• uj^'^'^f ^ therefore at zero 
temperature there is no dependence on A coming from G^^. We are interested in the poles 
of Gr at zero temperature. In this section UV constants will not be important. Theres 
are three characteristic regimes depending on the z/^^, the low energy [uj <^ fi) scaling 
dimension of the dual fermionic operator. The poles of the Green function are located in 
the lower half complex plane at 

u}c{k)=uj^{k)-iT{k). (3.19) 

For the three regimes we have the following [18], [22] 

• For (quasi-) Fermi liquid, z/^^ > i, 

io,{k) = vp{k - kp) + . . . , {k - kpf'^'F-^ ^ 0, (3.20) 

and the residue of the pole is Z = hivp. The pole represents a stable quasiparticle as 
one approaches the Fermi surface, with linear dispersion relation and vp being the Fermi 
velocity, vanishing decay width and a non-vanishing spectral weight Z at the Fermi surface. 

• For non- Fermi liquid, z/^^^ < ^, 

m 



ojJk) = (k- kp) '"''F , = const, (3.21) 

oj^{k) 



and the residue of the pole \s Z ^ {k — kp) "''f — > 0. The pole represents an unstable 
quasiparticle as one approaches the Fermi surface, with exponent in dispersion relation 
greater than one, the imaginary part is comparable to the real part of the pole, and a 
vanishing spectral weight Z at the Fermi surface. Non- Fermi liquid is example of a Fermi 
surface without sharp quasiparticle picks. 
• For marginal Fermi liquid, Uf^^ = ^, 

T,{(jj) ciu log m + idioj, ^ = —pi^, (3.22) 

where ci < and di are real constants. The single-particle scattering rate is linear in uj, 
while it is still suppressed compared to the real part as the Fermi surface is approached, but 
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the suppression is only logarithmic. The quasiparticle residue also vanishes logarithmically 
at the Fermi surface. We use this summary on quasiparticle poles below. 

Following [16, 17], we represent the fermion determinant in an effective action eq.(3.13) 
as a sum over poles of the retarded Green function in the black hole background. We obtain 
an analog to eq.(C.9) of one-loop action 

where z*[A(r)] is a functional of the order parameter A(r), V2 is the boundary spatial vol- 
ume. In order to make a connection with the field theory eq.(C.9), the following equations 
has been used for the complex frequency z^, 

\T{^+iz)f = — (3.24) 
'4 ^' cosh(7rz) ^ ^ 

It was shown in [17], that it holds for a complex z by matching poles and zeros of the two 
meromorphic functions. Here [F^ -|- ^-2)!^ = + iz)r{^ — iz). We will consider the zero 
temperature limit, therefore as explained in [16], the sum in eq.(3.23) is saturated by one 
pole eq.(3.19), 

z*[A(r)] = Lo,[A{r)] - iF[A(r)], (3.25) 

where the real and imaginary parts of the dispersion are functionals of the order parameter 
A(r). As was shown in [17], equation (3.23) for the fermion determinant captures only 
the singular contributions incorporated by the closest to a; = pole eq.(3.19), and smooth 
analytic terms are not improtant. 

We take the functional derivative, 

V2 ^f2Mr) , T\gH\ S , f I i-'. [A(r)] , 1 „ 



(3.26) 

with ^g^l^j.^ = 6{r — r') and therefore the dimension of the functional derivative is 
For the gap equation, = 0, we have 

AM - G^ntm^ ( 5u:.[A{r)] iz^Ajr)] 1 6T[/\{r)] iz.[A{r)] 1, 
^ > ~ TT TT V 5A(r) ^ 2ttT ^2' 5A(r) ^ 27rT +2- 

where ^ is the digamma function, '^[x) = ^^"^^^^^ • 
In the zero temperature limit, T ~ 0, we have 



(3.27) 



TT TT V ')'A(r) '2 a;^[A(r)]' 

^r[A(r)] 2vrT \ 

Generally, it is difficult to find the dependence for the pole z=i,[A(r)] [21]. However, here 
we have simplifications. First, the order parameter enters the Dirac equation essentially 
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as a mass term on the diagonal, i.e. it does not mix and ip as a. superconducting gap 
does. Therefore in the pole of the Green function eq.(4.17), the frequency is not affected, 
and there is only a shift in the Fermi momentum kp- Second, we can consider the order 
parameter A to be small. The procedure of finding the Fermi momentum is reduced to 
finding the bound state of the Schrodinger equation with zero energy [18] or to finding a 
solution to the Dirac equation which is normalizable at the boundary [14]. Since the order 
parameter contributes to the potential term in the Dirac equation, the shift in A;^? is given 
by the first order perturbative correction 



5kF[/^{r)] = / dr^V^0(r)taV°(OA(r), (3.29) 



kp kF±5kF[/^{r)] 
bhi 

where the zero modes tp^ are solutions of the (free, without A) Dirac equation with a; = 
and k = kp introduced in the next section, the signs it refer to components -F1/-F2 of 
the Dirac equation (A. 42), 6 is a dimensionless constant which has to be determined from 
the equation for the Fermi momentum kp in the presence of the gap A. Note that F* = 
diag(i(T^, icj^). The unperturbed Fermi momentum is kp = \J q^[i^ — and K is the AdS 
radius. In eqs.(3.20)-(3.22), the following substitution should be made k — \j2\qH\l (this 
substitution in the pole was shown to be true in [16] based on scaling arguments, see also 
Appendix A.l). Introducing magnetic field lowers the Fermi energy. For large magnetic 
field the lowest Landau level / = dominates and higher Landau levels are not important, 
while for small magnetic field all Landau levels should be included in the sum to correctly 
reproduce the limit of zero magnetic field. Following [16], we tune magnetic field to the 
point when k ^ kp, and once the pole crosses the Fermi surface it is counted in the fermion 
determinant of the effective action. 
For the Fermi liquid, u^p > \, 

a;4A(r)] = vp5kp[/^{r)l r[A(r)] ~ {5kp[/\{r)]f''^F . (3.30) 



Near the Fermi surface we have for small A ~ 

^a;4A(r)] .>i4..0(,)t^i..0(,)^ 6T\^ ^ f{r)}aH^ 



{5kp[/^{T)]) 



5^{r) i?3 ^ ' ' ^ ^ " 5A(r) R"^ 

(5/fcF[A(r)])2'^^F-i ^0. (3.31) 



r[A(r)] ^ ,,,_r.,„,n2.,,-l 



a;4A(r)] 

Therefore at T = 0, the gap equation gives the following solution 

A(r) = V iry'^ [r), (3.32) 

where Gint = -^t, and 6 is a dimensionless constant. The difference in factor 2 with the 
(2 + l)-dimensional case eq.(C.12) is due to taking one pole eq.(3.19) instead of two poles 
in the field theory, which does not affect our conclusions. Eq.(3.32) contains the radial 
profile of the order parameter, 'il;^[r)^a^ip^[r) shown in Fig. (2). The prefactor in eq.(3.32) 
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contains information about the magnetic catalysis for the fowest Landau level at T = 0. At 
zero temperature, the only solution is a nonzero gap which is proportional to the magnetic 
field and a radius of the four- Fermi interaction A ~ For a finite T, from eq.(3.27), 

there exist also a trivial solution A = 0, and a critical temperature Tc separates the phases 
with zero and nonzero gaps. This is in complete analogy with a (2 + 1) field theory case 
[1]. We consider the regime around Tc and the corresponding phase transition in the next 
section. 

For the non- Fermi liquid, v^p < ^, 



1 



u}4A{r)] ~ r[A(r)] ~ {dkF[A{r)])'-'^F . (3.33) 

Near the Fermi surface we have for small A ~ 
^^*[A(r)] .5r[A(r)] 

zeros for both term in the gap equation mean that there is no insability. Therefore in this 
case 

A(r) = 0, (3.35) 

i.e. no gap is generated in a magnetic field. Thus, for Fermi liquids, particle-hole pairing 
i/jTtp, is favorable in the magnetic field, while non-Fermi liquids do not support the pairing. 
The same conclusion has been reached for the case of the superconducting pairing < ^lJT■^p > 
where F contains F^ in [14]. There, it was suggested that taking the long-range four- Fermi 
interaction may generate the instability for non-Fermi liquids. On a technical ground, the 
momentum/frequency dependent four- Fermi interaction will change a simple shift in the 
Fermi momentum eq.(3.29), so that the derivatives of the real and/or imaginary parts of 
the pole will not vanish. 

4. Pairing instability in the Ginsburg-Landau formalism 

In this section we consider the pairing particle-hole instability using Ginsburg-Landau ap- 
proach in the bulk. Our calculation closely follows the leading order (one loop) Ginsburg- 
Landau procedure, with the only difference of using bulk fermion propagators. The pre- 
scription to construct propagators and vertices on the gravity side is given in [14], [22]. 

4.1 Microscopic calculations of magnetic catalysis. Ginsburg-Landau in a holo- 
graphic approach 

In the mean field approximation, the four-Fermi (contact) interaction ip{x)Tilj{x)il'{x)Tip(x) 
gives the following bilinear terms 

^j:T^pA + iJr^pA\ (4.1) 

where the second term is hermitian conjugate to the first one and F = F^FtF*, and the 
order parameter 

A = GinSTi;) (4.2) 
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is a singlet (number) in spin space. Here F = iV'^T^. The one- loop (Euclidean) action is 
given to the second order in A 

= [ [ d^xd^x'trg{x',x)rAix)g{x,x')rA{xy, (4.3) 

J ^(jint J 

where the Euclidean non-interacting Green function in the bulk is Q{x, x') = —{ip{x)'ip{x')). 
Assuming translational invariance along the spacetime {r, x}, we perform the Fourier trans- 
form 

g{x,x') = T^I ^a(r,r',ia;„,fc)e— (4.4) 

where the radial coordinate is r and the boundary directions are {r, x}; the fermionic 
Matsubara frequencies are = 7rT(2n-|- 1). We assume that the gap depends only on the 
radial direction, A = A(n). The one-loop effective action is given by 

5(2) + I dr'./g(?)A{r)A{r'rF{r,r')^ 

F{r, r') = -2TY,j (^2 *^^(^'' ^' ^)^'3{r, r', -ia;„, -^)f . (4.5) 

We make analytic continuation of the Euclidean Green function into the lower (upper) half 
plane in imaginary frequency plane, and use the following expressions relating Euclidean 
and retarded (advanced) Green functions [14] 

G{z)Q{-z) = g''{z)G^{-z) 

a^(r,r') = -a^(r,r')*, (4.6) 

where z = iujn, in order to rewrite the action in terms of retarded (advanced) Green 
functions 

(2^ y — tanh — trg«(r', r, A;)*Fg^(r, r', -Q, -k)T, (4.7) 

where we substituted the Matsubara sum by the contour integral, iun — )• z, and Q is on 
the real axis of z. In order to calculate this integral, we express the bulk Green function 
through the boundary one as given in [14]. The bulk Green function is a solution of the 
free Dirac equation, eq.(A.32), 

D{n, k)g^{r, r', k) = -^i6{r, r'), (4.8) 

with the free radial Dirac operator D{Q, k), which includes the mass term and the magnetic 
field but has zero gap, A = 0. The bulk Green function is constructed through the modes 
■0(r), ip = Q-^^i+^kx^^^^^ which are solutions of the free Dirac equation eq.(A.32) 

Din,k)^Prad^al{r)=0■ (4.9) 
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Due to the choice of the Gamma matrices, eq.(A.15), V decouples into two-component 
spinors, ij) = ('i/'i, ^^2)^) which are eigenfunctions with definite eigenvalue of r'T*r^. There- 
fore the bulk retarded Green function has the block-diagonal form, eq.(A.40), where the 
components Qa, ol = 1,2, are constructed from the solutions of the Dirac equation as [14] 

n ^, _ G^{n,k) i (r)^»(rO r>r' 

Qa {r,r)- — ^ X I ^ ^ ^, , (4.10) 

with 'ipa = iip'^a^. Note that F'^ = F* = diag(i(7^, ia^) and the dimension of ■0 is ["0] ~ L^'''^; 
the minus sign comes from the definition of the bulk Green fuction Q = — < iptjj >. Here 
the prefactor arises from the Wronskian [14], and includes the retarded Green function of 
the boundary field theory (later we refer to it as the boundary Green function). Since the 
Wronskian is a constant related to the conserved charge current, it is simpler to calculate 
it at the conformal boundary [14] with the result given by eq.(4.10). 

There is the following reasoning behind contruction of the bulk Green function in 
eq.(4.10). At the boundary, r,r' — )• 00, the behavior is given by two terms (we omit 
prefactors) 

G~r^'^-^ + r-^'^, (4.11) 

and we put to zero the non-normalizable part ~ qt^^ leave the normalizable part 

~ r~^*. At the horizon, r,r' — >• tq, there are two terms 

G~ e-^*^''^ +e^'=^o, (4.12) 

where we throw away the outgoing solution e~*'^''°, and leave the ingoing one e*'^^". In this 
way we obtain the retarded Green function. Note that e^*'^''" ~ (r — tq)^'^^/'^ with the 
Lorentzian signature. In the Euclidean space, this means that when the behavoir of the 
Green function is fixed at the horizon, r, r' — )• ro as G ~ (r — r^)^^/'^ , then the retarded 
Green function has asymptotics Gr ^ {r — ro)^/'^ , with u > (and corresponds to the 
ingoing solution), and the advanced Green function has asymptotics Ga ~ (r — ro)~^/'^ , 
with a; < (and corresponds to the outgoing solution). 

The boundary conditions are fixed as follows for the solutions of the Dirac equation. 
The (free) solution of the Dirac equation has the following behavior near the boundary, 
r 00, eq.(A.39) with A = 0, 

i^a ~ a«r™^ 1^ ° ^ + 6„r-™^ |^ J ^ . (4.13) 

The two spinors (1, 0) and (0, 1) are eigenstates of with opposite eigenvalues, implying 
that and ba are canonically conjugate (in a radial Hamiltonian slicing) [23]. Therefore 
a boundary condition must be imposed on one, with the other allowed to fluctuate. For 
mR > I , we choose the fluctuating piece to be the normalizable mode at the boundary (with 
regular behavior) proportonal to (1,0). This gives us the solution tj;^'^^. More generally, 
the quantization choice for mR > | is to impose the boundary condition = on the 
fluctuating mode. Another quantization choice for mR > 5 is discussed in [23, 18]. For the 
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solution V'*" we impose ingoing boundary conditions at the horizon, r = tq. Thus, there are 
two normalizable solutions with the following behavior at the conformal boundary r — )■ cx) 




C" = ^r-«( : )+r-"^ ( ; ), (4.14) 

where the boundary Green function Ga was introduced in ip"^^ following its definition be- 
ing proportinal to hajcLa in the ingoing solution. Using representation eq.(4.10), we can 
show that g^{r,r',n,k)* = -^^(r, r', J^, fc), since il}'"^y * = ip'^^y and V'*" * = V'""*- Also 
g^{r,r',n,k) = r*g^(r',r,17,fc)tr*, since ia^G^^ir' ,r,n,k)Ucr^ = -G^{r,r' ,n,k)* . There- 
fore we can rewrite the kernel eq.(4.7) in equivalent form 

7^ / — tanh— trr*g^(r,r',J7,A:)tr*rg^(r,r',-0,-fc)f. (4.15) 

Using the relation eq.(4.10), we obtain 

72 u roo 



^^"'"') = i^ j W? j J^'^''^^G,{ii,krG,{-Q,-k)x 



{ i)f{r\ n, k)^(j^%l^f{r', -n, -k)ijfy{r, f], A?)ttTVf ^(r, -k), r > r' 
\ ^l^^{r', n, fe)t(7V2'^^(r', -S^, -k)ipY'{r, n, fc)tc7Vr(r-, -k), r < r' 
+(1^2), (4.16) 

where we used V'q"!"^) = ^a"* '0a"'(^)* = which follows from the definition 

of in(out)going solution, the solution ^'"^'^ is real, and il^i{—k) = il^2{k) which follows from 
the symmetry of the Dirac equation [20]; we omitted other arguments by ip^s. 

All quantities in eq.(4.16) can be obtained numerical. However, at low temperatures 
T <ti fJ-, calculations can be done analytically, due to the fact that the main contribution 
comes from the pole in a retarded Green function describing the Fermi surface, i.e., the 
closest to the origin pole ~ as ~ fei?. The same argument that physics occurs close 
to the Fermi surface in the bulk was used in [16] describing phenomena in magnetic field 
and for the BCS theory in a holographic approach by [14]. This repeats the reasoning of 
the field theory, as in the BCS theory. Close to the Fermi surface and at low temperatures, 
i.e. at T <^ 0, <^ 1^1, the boundary Green function, G = Gi, to the leading order is given 
by [18] ^ 

where k± = k — kp is the perpendicular distance of the momentum from the Fermi surface, 
hi and Vf are real constants and are calculated below, /12 is positive and the phase 9 is such 
that poles of the Green function are in the lower half complex frequncy plane, i' is the zero 
temperature conformal dimension at the Fermi momentum, i' = Vkp^ given by eq.(A.46) 
with A = 0. 
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We perform the momentum and frequency integrals in eq.(4.16). We make the same 
assumption as in [16] and as in our variational calculations. We consider only contributions 
near the Fermi momentum k ^ kp (see discussion after eq.(3.29)). Therefore, in the 
leading order, there is no momentum dependence in the boundary Green functions, and 
the momentum integral (Pk is trivially performed, which gives in the magnetic field a factor 
\q%\. In the frequency integral, for T ^ ri, we substitute tanh — 1, and we have 

1 ^ hi hi 



, 1 /-"^ d\ 



2 -n/vF + h2e-^^+''^''n'^'' n/vF + h2e'^+^'"'n'^'' 

\ iljf\r', ^)^a^^ljfy{r', -Q)ijf{r, Jl)t(jVf (r, -Vi) r < r' ' ^ ' ' 

where the wave functions are evaluated at the Fermi momentum. In the frequency integral, 
depending on the critical exponent z/ either the first or the second term in each denominator 
dominates. To make this comparison, note that /12 has a dimension, i.e. /i2 ~ fi^^'^'^ . At 
small frequencies O ^ ;U, for v > ^ the first term ^ Q. dominates, while for u < ^ the 

second term ~ ( j dominates, as also shown in [18]. As discussed in [14], the range 
T <^ 0, <^ fi implies that the fermion wavefunctions should be evaluated at = and 
k = kp in the extremal T = black hole background. This means, they are exactly zero 
modes at the Fermi surface, i/^^'^^ = = which will be calculated below. 
For u > the leading behavior of eq.(4.18) at small temperatures T <C /U is 

F{ry) = -^H?^ V°(r)VV°(r)V^°(rOtaVO(r'), (4.19) 

where ijj^ is the T = fermion zero mode at the fermi surface. 

For u < ^, the leading behavior of eq.(4.18) at small temperatures T ^ /_f is 

^^"'"'^ = "^3i^'°'^^^^^^^~T^^ V'^^taV^OV'^rOtaV^r'), (4.20) 
that for I' < \ becomes to the leading order in expansion T//.i temperature independent 
Fir,r') = -|^|L^cos(2vrz.)^^ V'^^taV^OV'^rOVV^r'), (4.21) 

and for \ < < ^ has a wrong sign for the kernel F to give a nontrivial solution of the gap 
equation. This means that there is no instability in particle-antiparticle pairing for u < ^. 
We observed this already before in variational calculations, where non-Fermi liquids did 
not support ^'''/'-pairing while the Fermi liquids did, see eqs.(3.32) and (3.35). 

For the Fermi liquids, we obtain the equation for the critical temperature from the 
effective action eq.(4.5) 

^ + / dr'^/^^)A{r')F{r', r) = 0, (4.22) 
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which is an analog of the gap equation but with free propagators (Green functions). As in 
[14], we use the factorisation of F{r,r') to write an ansatz for the gap function 



A°(r)~ V°(0"^t^V°(0- (4.23) 
Using this ansatz back in the equation, we get the following critical temperature for f > \ 

Tc = ^""l^'.^e^'"^ / drv^(^°(r)taV°(r))^ (4.24) 

The radial form of the gap function eq.(4.23) was obtained by variational calculations 
in eq.(3.32). Eq.(3.32) can be considered as a one loop mass gap equation (self-energy 
correction) with a four-Fermi vertex Gint and the bulk fermion propagator Q^{r, r') given 
by eq.(4.10) and modified as in eq.(3.29) to include the gap. The mass gap equation reads 
A(r) ~ Gint J (-2^p / d(jjG{oJ,k,r,r), where integration gives ~ \q7i\hiVF and the resulting 
gap as in eq.(3.32). Substituting eq.(3.32) in eq.(3.29) for Skp, we have 



GintVp 

IF 



^eff ~ ^-^^dr,/^){tl^\r)^a'^p'ir)f, (4.25) 

where N^jf is the effective density of states at the Fermi surface. 

In the next section we calculate the constants hi and vp for the different fermion 
charges q, and find the zero mode wave function tl;^{r). We follow procedure outlined 
in [14]. We then find the critical temperature as a function of the charge q and the 
magnetic field. 

4.2 Solving for the zero modes and finding the critical temperature 

In this section we follow the procedure of [14]. We solve the Dirac equation for the zero 
mode in black hole background, that amounts to finding a solution at zero frequencies 
a; = (we use here cj = fi) in the T = background. We will also set mass m = 0. As 
was shown in [14], analytic solution can be found in this case. The Dirac equation in the 
magnetic field 

Vm^ = 0, (4.26) 

where Vm = du + \^abM^"'^ — H^M is defined in Appendix A.l. The Dirac equation can 
be written explicitly including the spin connection as 



dr + Jgutcr in ^=0" (w + qAt) H — ^=0" - 

^ \/-gtt V-9tt 2 




- cr^^^xfx - (^^^^yfy - A ) 1 ( ) = 0, (4.27) 

where A — ?• ^^/2\q7^\l is the eigenvalue for the Landau levels which takes into account {x, y}- 
plane physics in the magnetic field (Appendix A), and tf^ = (■0i,V'2)"^; and the F matrices 
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are given by eq.(A.15). In the basis eq.(A.15) the two components decouple, therefore 
below we solve for the first component. Substituting the spin connection we have 

~^^dr + - -^cT3(a; + qA,) - ^^^(3 + ^) - a) F = 0, (4.28) 

with F = {yi,y2)- As in the AdS2, it is convenient to change the basis eq.(B.14) 



y 1 \ / 1 -i\ ( yi 

m \ -i I \ y2 



(4.29) 



that simplifies the second order differential equation for one component. The Dirac equation 
is given 

'"^-a'dr + ^ia^m + -^a\oo + qAt) - a^!^(3 +'^)-x]f = 0, (4.30) 



i?2 R ' ^ \ ' 2i?2 ^ 2/ 

with F = {yi^my. 

We introduce dimensionless variables with the goal to scale away the ^^5*4 radius R 
and the horizon radius ro 

m 
K 

M^r^M, Q^rlQ, ,H^rlH (4.31) 

and 

{t,x) ^ —[t,x), Am ^ ^Am, i^^-^oj, 

ds^ R^ds^. (4.32) 
In the new variables we have 



At = - ^), 11 = V^gprl (4.33) 
r 



and the metric is given by 



^ r\-fde + dx') + \^, (4.34) 



with the horizon at r = 1, and the conformal boundary at r — )• oo; the red shift factor is 

The Dirac equation is given by 

(^-r^^/fcT'^r + ria'm + ^^'(^ + 'Z^*) " ^'^(^ + ^) - F = 0. (4.36) 
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We set m = 0. Then we get the following second order equations for each component 

r'fd', + {f>r'f + r'f')dr + ^-^r'f + 2r'f' + ^ 

+ ^((a; + qAt) ± '-^f ^ r^iqA', - A^) y,.^, = 0, (4.37) 

with At = /x(l — ^), and the upper/lower sign is for yi/y2- At T = 0, from eq.(4.33) we 
have r*^, = 1, and the red shift factor develops the double zero near the horizon, 

f=ir-mr^^^-^2r + 3)^ ^^^^^ 

Due to this fact, the metric near horizon reduces to AdS2 x R^, and calculations in AdS^ 
are possible to do analytically at small frequencies [18], [14]. We will utilize that below. 
We introduce a new radial variable z, 

1 



1-z' 

then the second order differential equation is given by 



(4.39) 



+ + qAt) ± ^)2 ^ iqA't - A^j y,,2 = 0, (4.40) 



with 



f = 3z^{z-zo){z-zo), zo = ^{4 + iV2), 
At = fiz, 11 = VSgprl, (4.41) 

and horizon is at z = 0, the conformal boundary is at z = 1. Due to the double zero of 
the red shift factor near the horizon, the second order differential equation eq.(4.40) can 
be solved analytically at oj = [14]. 

We put for completness the second order equations for the components when spin 
connection has been eliminated using transformation eq.(A.8) 

= {r^Vfy'^^'^, (4.42) 

where equations for ^ contain the spin connection, and for $ do not. The second order 
equations in dimensionless variables without spin connection read 

(rV9.^ + (2r3/ + + !^ + i(a; + ± ^ r\q^t - A^) yi:2 = 0, 

(4.43) 

and we use the same notations for the components for $ = (^1,^2)"^ as for ip] again the 
upper/lower sign is for yi/y2- Using the radial coordinate z, these equations are written as 

+ ^d. + {^ + j{u; + qAt± ^)2 T iqA't - A^) y,,, = 0. (4.44) 
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These equations look simpler than the corresponding equations with spin connection. How- 
ever, they are not readily recognized by MAPLE program, and we will not use them. 
Writing the transformation which removes the spin connection explicitly as 

-1/2 



^ = (^ (i_^)3 V3(z-zo)(^-^o)j • (4.45) 

gives an idea about the prefactor that we should expect to get in the solution. As men- 
tioned, these equations are not recognized my MAPLE. We therefore proceed with eq.(4.40) 
which contains spin connection. 

Near the horizon, z = 0, we have / = Qz'^ and 

Gz'y" + I2zy' + + M! _ a^)^ = o (4.46) 
and the same for z, giving the behavior near horizon 



5±^ ^ = ^V6A2-(g^)2, (4.47) 



2/1 ~ 2/2 ~ ^ , - „ 

D 

with /Li = V^gpi"*- We will be interested in this scaling exponent given at the Fermi 
momentum, 

i^^iykp = \^Qkl-^q^glrt (4.48) 

which is the conformal dimension found in eq.(A.45). Note that r^^ = 1 at T = 0. 

Putting (x) = eqs. (4.40,4.41) and using MAPLE, we find the analytic solution for 
the zero mode [14] (see also [24]). The first solution with regular behavior z~2+'^ at the 
horizon, z ~ 0, is given by 



i(-lTV2W^o) 

y% = Ni.^{z-1)-2z-^^^^^{z-zq) 




z - zo 

, . ^/2 qa 2iV2z \ 

, I _ + =f ± i|, 1 + 2n, ^^^^^ j , (4.49) 

where 2F1 is the hypergeometric function, Ni,N2 are normalizations defined later, and 
upper /lower sign is for 2/1/2/2; the role of momentum is played by A — t- y^2\q'Hl\. The 
second solution, with behavior z~2~^ at the horizon, is obtained by replacing vx — )• —i/x 
in eq.(4.49) 

€,2 = Ni;2 (j^^ with ux ^ -u^ , (4.50) 

and it will be required to have a regular behavior at z ~ for small frequencies. Since 
normalization factors are constants, we find their relative weight by substituting solutions 
back into first order differential equations eq. (4.40,4.41) at z ~ 0, 

A^i 6iux + q^l fzoY^^'""' iVi ^ Gzi^A^ /^y"/^^° ^^^^^ 



N2 VqX \zoJ ' N2 \/6A \zoJ 
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Figure 1: Fermi momentum kp vs. charge of the fermion field q' = \/?>q. We choose = 1, 
= 1, therefore /i — \/3- The inner (closer to x-axis) line of the filled/shaded area is Vkp = 
and the outer line is Vkp = 5 1 so that the shaded region corresponds to < Vkp < 5 • At a given q 
there are multiple Fermi surfaces. From left to right are the first, second etc. Fermi surfaces. They 
dissapear at v^p = 0. Positive and negative kp correspond to Fermi surfaces in Gi and G2 Green 
functions respectively. At g = 1, which on the plot is g'-\/3 kp ~ 0.92 in agreement with [20]. The 
first Fermi surface hits the border-line between a Fermi and non- Fermi liquids v — ^ &t q' — 2.71. 



The same relations are obtained when calculations are done for any z. The zero mode 
equals to V'" = \{yi + im, m + m) with §1-2 = §1-2 + ^i;2- 

To obtain the Fermi momentum, we can follow [18] and search for a normalizable 
solution of the Dirac equation with the certain boundary conditions, which is equivalent 
of looking for a bound state in a zero-energy Schrodinger equation. We use an alternative 
way suggested in [14], that also uses that the solution should be regular at the horizon and 
obey certain falloff conditions near the boundary of AdS^. To construct Tp^'^y, we required 
the boundary condition a = on the fluctuating mode, i.e. at the conformal boundary 
z 1 

^ \y2 + wi ) \ y 

Therefore the equation for the Fermi momentum kp \s 

\i-miz-l)-'^/'^{y2 + iyi) = Q. (4.53) 
Using the zero mode solution eq.(4.49) in the eq.(4.53), we have 

2Fi(l + + '-f)^2+ ^fcF - 1 + 2z^fcF> 1(1 - i^) _ <^^k^ - nk^ 



2Fi(z.,, + M^, 1 + _ v|?M, 1 + 2vu,. lix _ ,V2)) /cF(-2^ + ^2) ' 



(4.54) 



with V].p = 6A;^ — isil^)^- We solve the equation for the Fermi surface numerically, using 
MATHEMATICA to evaluate the hypergeometric functions. The solutions of eq.(4.54) are 
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depicted in Fig.(l). There are multiple Fermi surfaces for a given q. Following [14], the 
largest is called the first Fermi surface, the next \kp\ the second Fermi surface, and so 
on. For all further plots, we choose r^, = 1 and gp = therefore ^ = \/3. We recover a 
result of the numerical solution of the Dirac equation [20]: for q = 1 which is = ^/?> we 
have kp = 0.9185. In Fig.(l), positive and negative kp correspond to the Fermi surfaces 
in the Green functions Gi and G2. The relation between two components when m = is 
— Gi(L,fc) I^*-*]' therefore Fig.(l) is not symmetric with respect io kp — Q axis. 
We substitute the Fermi momentum into the zero mode solution eq.(4.49) and get 
the radial profile for the pairing gap function A" given by eq.(4.23). We plot A'^(z) for 
different charges. Fig. (2). The curves are normalized to have the same maxima. Charges 
are increased from left to right. For large charge, when 1/ > ^, the zero modes are supported 
away from the horizon, while at smaller charge, when i/ — )• |, the zero mode functions are 
supported near the horizon. The same tendency was first observed for the Cooper pairing 
[14, 21]. This means that for non- Fermi liquids, at u < ^, the physics of the Fermi surface 
is captured by the near horizon AdS2 x region. 

gap 




z 



0.2 0.4 0.6 0.8 1.0 



Figure 2: Wavcfunction of a pairing mode A° — ip'-^'^ a^ip'^ as a function of the radial coordinate z, 
with the horizon at z = and the boundary at z = 1, for different values of the charge q' = \/3q 
for the first Fermi surface. We set ^ 1, gp ~ 1. From left to right the values of the charge are 
q' = {3,3.4,4,6,8,10}. The curves are normalized to have their maxima all the same. At small 
charge, non-Fermi liquid, the wave function is supported near the horizon. At large charge, Fermi 
liquid, the wave function is supported away from the horizon. 



We express the boundary Green function through the zero mode solutions. To do 
that we should obtain the solution at small but nonzero frequency and expand the Green 
function to the leading order in cj. We follow the matching procedure of [14] between the 
solution in the "near" (to horizon) region ("inner" in terminology of [18]), z <^ 1, and in 
the "far" (from horizon) region ("outer" in terminology of [18]), z ^ oj. Matching the two 
solutions gives a solution on the full spacetime as long as the near and far regions overlap, 
a; <Cl. 
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Near the horizon, z <^ 1 f{z) = Qz^, and the second order wave equation, eq.(4.40), 
becomes 

z^a^ + 2zd, + 4 + -f ± 6i- + 9 - 6X^)] yi;2 = 0. (4.55) 



36 v ■ z' z 
Using MATHEMATICA, we obtain the following solutions 

= Ci;2Z-i--^e-^iFi Q ^ 1 + i.;, + ^, 1 + + Z)i;2(^A ^ -i^A), (4.56) 

with upper/lower sign is for yi/y2- Requiring that the solution is ingoing at the horizon 
2 ~ 0, ~ e+*'^/6^ fixes the ratio 



Ci C2 T{-2vx)T{l + vx-'-f), . .2u, 

— — [ — lUJ)^. 



(4.57) 



The near horizon solution is in the matching region z ^ oj, 

= A,,2Z-'^-''' + B,,2Z-'2+'^\ (4.58) 
that corresponds to the AdS2 boundary, given in eq.(B.24). From eqs.(4.58,4.57) we have 



G^^^^^^^ r(-2.A)r(i + ..-f^) 



(4.59) 



Ai A2 r(2i.A)r(i-i/A-T) 

since by definition the ratio Bi/Ai gives the retarded Green function in the IR CFT living 
on the boundary of AdS2. We calculated it in the Appendix C, eq.(B.30). 
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Figure 3: Constant hi, reflecting the UV physics of the AdS4 bulk, vs. charge q' — ^/3q. It 
vanishes at Vkp = 0. The multiple lines are for various Fermi surfaces, in ascending order with the 
first fcrmi surface on the left. Note, hi has the same sign as kp- As above, positive and negative 
kp correspond to Fcrmi surfaces in the Green functions Gi and G2 respectively. 



In the asymptotic far region z ^ uj, the second order wave equation is solved pertur- 
batively in to, 

y(^ = yn+uy^l (4.60) 
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where ^(^^2 includes the zero modes found m before, yf'l = yi.2 + ??i.2- Expandmg 

y(0) in 

the matching region, z <^ 1, 

S,.^2{l^) = (-l)3/2(_^„)-|+.. (^fO j ' . (4.61) 



Comparing the near solution eq.(4.58) and the far solution eq.(4.61) in the matching region 
to <^ z <^ 1, we get 

{-zoY^^^G^^, (4.62) 



that determines the relative contribution of and eta". We also have relations for 
and given in terms of A'^i and A'^i, eq.(4.51). Comparing eq.(4.59) and eq.(4.62), a; and 
z scale with the same power in the solution y around the horizon. From eq.(4.62) follows 
that 

r??.2 ~ ~ G^-^ ~ w^^^ (4.63) 

The first order correction y^)^ satisfies an inhomogeneous second order wave equation with 

y^\ as the source. To calculate the retarded boundary Green function, we need only the 
leading asymptotic behavior near the boundary z — t- 1. The asymptotic behavior can be 
found by integrating the Dirac equation (4.40) as in Appendix C of [18] with the result [14] 

g^" + «i" = 2i(l-.)» ^°''''^<l£+fe°l^' . (4.64) 

Vi W2 

Note, that from eq.(4.64) and eq.(4.63), only fj depends on the frequency, while all other 
wavefunctions are independent of u. 

The Green function of the dual field theory defined on the boundary of AdS4, as 
introduced in eq.(4.14), is 

G = lim^^±^. (4.65) 

m + lyi 

Expanding in small frequency uj, 

G = lim yUw| + vU^fi'2 (^gg) 

+ ij/O + 77° + i7?o + uj{y^2 + ^^i ) + 0{uj'^) 

where the zero mode wavefunctions yi.2, Vi-2 defined in eq.(B.24) with normalization 
given in eq. (4.51,4.62) and asymptotic behavoir of the last term in denominator given in 
eq.(4.64). Near the Fermi surface, k± = k — kp, and at T = the Green function is written 
[18] 

G = . 2. . (4.67) 

which was used in our calculations for Tc, eq.(4.17). Note that all quantities here In the 
above formula, the last term in the denominator comes from ffl.2 and includes G^^ ~ u^^''f , 
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Figure 4: Constant vp, reflecting the UV physics of the AdS^ bulk, vs. charge q' = VSq. It 
vanishes at — \- For the first Fermi surface it happens [vkp — \) q' = 2.71. The multiple 
lines are for various Fermi surfaces, in ascending order with the first fermi surface on the left. Note, 
vp has the same sign as kp. As above, positive and negative kp correspond to Fermi surfaces in 
the Green functions Gi and G2, respectively. 



i.e. it is determined by the IR AdS2 physics near the horizon. Other terms are determined 
by the UV phyics of the AdS4 bulk. The constants hi and vp are 

hi = lhn y'^'^' , (4.68) 



VF 



l(^'..v^<./)"i..ML±!|!, 

where all wavefunctions are evaluated at k = kp; the zero mode has components ip^ = 
{yi + iy2, 2/2 + ^y?); ^iid the wavefunctions yi.2 are given by analytic expressions, eq.(4.49). 
The constants /i2, can be also obtained from this expression. The constants hi and vp are 
dimensionless in eqs.(4.68),(4.69), i.e. the scaling is (from dimensional to dimensionless) 
hi —7- ^^1, Vp — 7- ^vp. Both constants are real. They are plotted as function of qfi in 
Figs. (3), (4). The Fermi velocity vanishes at the horizontal line vp = when z^^^ = |. The 
wavefunction renormalization hi vanishes when u^p = 0. The multiple lines in each plot are 
for various Fermi surfaces, starting with the first FS at the most left. Positive and negative 
Vp, hi correspond to the Fermi surfaces in the Grenn functions Gi and G2 respectively. 
Both Vp and hi have the same sign as kp. The Fermi velocity decreases at small charges 
and tends to the speed of light at large charges. Geometrically this means, that as the 
charge is lowered the zero mode wavefunction is supported near the black hole horizon, 
Fig. (2), where the gravitational redshift reduces the local speed of light as compared to the 
boundary value. This was observed also in [14, 18]. 
The dimensionless critial temperatute is given by 

r. = / d.V=^(V'°(.)taV°(^))^ (4.70) 
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Figure 5: Critical temperature vs. the charge q' — ^/iq for the first Fermi sm^face. Only 
parametric dependence is shown, and the I/tt^ factor is not included. We set — l^gp — 1. Note 
that Tc vanishes around q' ~ 2.71 which corresponds to v^p = \ for the first Fermi surface. This 
plot illustrates that pairing is supported only for v^p > | which is the region of Fermi liquids. 



with transformation to the dimensionless variables given by 

(2 \ ~2 2 
Gint, n^^n, hi^^hi, ^o^i?3/2^o. (4.71) 

Using the analytic expression for the zero mode eq.(4.49), and the results for hi^vp^kp 
calculated for a given charge q, we plot Tc as function of charge for the first Fermi surface, 
Fig. (5). The next Fermi surfaces give smaller contributions and are not depicted on the 
plot. The critial temperature vanishes exactly for v = which is for the first Fermi surface 
is at g/i = "v/Sq = 2.71 and with kp = 1.65. Therefore there is no pairing for v <\. This 
happens due to the fact that to the leading order the density of states eq.(4.25) vanishes 
for the non-Fermi liquids. This conclusion agrees with our variational calculations, where 
pairing occurs only for the Fermi liquids, while non-Fermi liquids do not support pairing. 
Note, that in principle fermions from different Fermi surfaces can participate in pairing. 

Now we analyze behavior of the boundary field theory with magnetic field. The con- 
formal dimension of the fermionic operator O in the IR CFT eq.(A.46) is 



where it is taken in the chiral limit m = and at the Fermi surface kp = q^i, and the 
dimensionless variables introduced before are used. We rewrite this expression in terms of 
Q and H, using eq.(2.8), as 




, kj, A2 g2g2 

Vl + (H/Q)' ^ 6 12{1 + {H/Qn ^■^^> 
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Figure 6: Conformal dimension of the fermion operator in the IR CFT Vkp as a function of the 
ratio ^ for different values for q' . We used parametrization v ^ \J a/ vT+ x"^ + h — c/ {1 + x"^) with 
X = ^. The curves from right to left correspond to increasing charges g', a — {0.3, 0.5, 0.8, 1, 1.1}, 
while h and c = 1 are kept fixed. Top plot is for 6 = 0, zero gap, and bottom plot is for h = 0.5, 
nonzero gap. At small and intermediate charges, the Fermi liquid regime and hence the particle-hole 
pairing are realized at a threashold value of the magnetic field, which is consistent with magnetic 
catalysis in graphene with impurities. 



Writing Vkp as a function of x = H/Q, v ~ y + ^ ~ i+x^ ; P^o* Vkp vs. x for 

different parameters a,b,c, Fig. (6). Multiple curves correspond to increasing values for a 
starting from the right to left curves, while b and c are kept constant. This corresponds 
to increase of the fermion charge q. At small charges the Fermi momentum kp is small 
(see Fig.(l)), that corresponds to small a (curves to the right). All curves show the rapid 
growth at the beginning, then saturation at some maximum, and fall off which happens at 
large enough x where we are not interested any more. For small charges, the rise to = ^ 
happens relatively quickly, and means there is a narrow window of magnetic field for non- 
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Fermi liquids. However, the region for ^ < u < 1, the case of Fermi Uquids, corresponds to 
a much wider range of the magnetic field. For small and intermediate charges, regime of the 
Fermi liquid requires a threashold magnetic field, when = ^ is reached. This is consistent 
with the magnetic catalysis in a system with dissipations, e.g., graphene with impurities, 
where there is a nonzero width and the particle-hole gap is induced at a threashold value of 
the magnetic field [1]. For large charges, corresponding to curves to the left, there is only a 
regime of Fermi liquids (no non- Fermi liquids are possible) and pairing is supported for any 
magnetic field. Top and bottom pannels are plots for zero and nonzero gaps, respectively. 

In this paper we have imagined the existance of an "experiemntal knob" which can be 
used to adjust the UV scaling dimension of fermionic operator. As noted in [22], the most 
useful knob will likely depend on the UV geometry into which this AdS2 is embedded. In 
our case, an external magnetic field Ti will allow one to tune the IR scaling dimension and 
to explore different sectors of the boundary field theory. 

In analogy to the superconducting instability < ipip > of the black hole, we find that 
the Breitenlohner-Freedman bound can also be broken in case of < ipip > condensate for 
large enough fermion charges q. Using the dimensionless variables, the conformal dimension 
of the bosonic operator in the IR CFT in the presence of a magnetic field is 

(4.74) 

which is obtained from eq.(A.45) by k"^ — )• 2\qH\l, the mass gap A — )• m with ni being 
the mass of the bosonic field living in the bulk. Here the last term ^ distingushes the 
bosonic case from the fermionic one. The conformal dimension of the bosonic operator in 
the UV CFT, A,^, in the dimensionless variables is given by 





m2 + . (4.75) 

Breitenlochner-Freedman (BF) bound is broken when conformal dimension becomes imag- 
inary. As in the case with superconductor, there is a parameter range where the conformal 
dimension in the IR (coming from the near AdS2 horizon) is imaginary and the conformal 
dimension in the UV (coming from the AdS/j^ bulk) is real. Expressing the mass m through 
A^, the condition breaking the BF bound is 

qVri- > 2A0(A^ - 3) + 3 + ^^Mi. (4.76) 



Here we restored the dimension. 



5. Equation of state and transport properties of the boundary field theory 
at zero magnetic field 

In this section we consider thermodynamics of the boundary field theory, namely we obtain 
an equation of state and find the scaling behavior of the specific heat with temperature. 



-28- 



Then we consider transport properties of the system on the boundary, specifically we 
calculate the DC conductivity and analyze its scaling behavior. We do not specify the 
boundary theory. Instead we use the "dressed" by the gravity fermion propagators obtained 
from the AdS3+i/CFT2+i analyzes in [18]. As a result we obtain behavior of systems with 
properties ranging from Fermi and marginal liquids to non-Fermi liquids. In particular we 
reproduce correct temperature scaling for the DC conductivity and specific heat in case 
of the Fermi and marginal liquids. Since the two-point Green function of the boundary 
theory has been obtained using the AdS/CFT correspondence, it is "exact" in terms of 
gauge coupling corrections. Therefore the lowest order diagrams on the field theory side 
should suffice. Of course we lack the knowledge of the "dressed" by the gravity gauge- 
fermion vertex. For the quantities considered below, the scaling behavior does not change 
when vertex corrections are added. 

5.1 Equation of state and specific heat 

An effective potential in the CJT formalism is given by [25] 

r,ff = ^Trln5"i + ^TriS^'S - 1) + T2[S], (5.1) 

where S* is a dressed fermion propagator, r2 is the sum of all two-particle irreducible (2PI) 
diagrams, and trace Tr involves integration J d?x. The last two terms can be simplified 
with the help of Dyson-Schwinger equation, to give 

re// = ^Trln5-i- jTr(S5), (5.2) 

where the self-energy is E = 5"""^ — S^^ . 

We use "dressed" by the gravity retarded and advanced fermion propagators [18] 

n ( Z\ i-hivp) 



uj — VFk± + S(a;, kp) ' 

GAioj,k) = -GRiu,kr = i'^iw ^ V (^-3) 

CO — vpkx + E*(a;, kp) 

where the momentum is counted from the Fermi surface k± = k — kp, hi and vp are 
real constants (we keep the same notations for the constants as in [18]). The self energy 
E = hvFc{kp)oj^'^*'F contains the real and imaginary parts, E = Ei -|- zE2, with imaginary 
part coming from scattering processes of a fermion in the bulk, e.g. via pair creation, and 
scattering into the black hole. The spectral function defined as A(oj,k) = ^\m.GR{uj,k) is 

M-rk) = --, — 'r,^^^'^^ , , , (5.4) 

t: {uj - vpk^ + T.i{oj,kp)f + i:2{oj,kpY 

Exactly due to inelastic/dissipative processes we are able to calculate transport coefficients, 
which will be infinite otherwise. However, the imaginary self energy gives rise to a branch 
cut in the fermion propagator along Imw = in a complex u plane [16, 26, 27]. Therefore 
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in calculation of Matsubara sum we should take into account contributions from poles and 
discontinuities along branch cuts [26, 27] 

T Y Pi^^rn) = Y n{zi)Res{F, z = z^i)-Y T ^MC)Disc F, (5.5) 

oddm poles cuts'' 

with analytical continuation iujm — )• -z, and n{x) is the Fermi distribution function. One can 
use either n{x) = n{^) or tanh(^) functions with prefactors (—2^) and (—4^) respec- 
tively in the contour integral which give the same result for the observables. Calculation 
of Matsubara sums using perturbative expansion in the imaginary part of the self-energy 
has been developed in [28]. 

For simplicity we take {—hivp) 1 which will not change our results qualitatively. 
Using the retarded fermion propagator eq.(5.3), an effective potential is 

I V2 f (Pk f ^ z 
^eff — ^ — — -7^ / 7 — TFT / dz tanh — - X 
Ani T J (27r)2 2T 

rj^fl. z - VFkj_ + T,{z,kF) 1 T.{z,kF) \ , . 

^2 T 4z-VFk^ + J:iz,kF))' ^^-^^ 

where we substituted Matsubara sum by the contour integral. The original contour Co 
going around the poles along imaginary z-axis was deformed in the contour C going along 
the real z axis and F being arcs at infinity with vanishing contribution [26]. In case of a 
real self-energy the result for the contour integration is (see Appendix C and [16]) 

I (i^E(5^>"(2-h|^) + ^E(..)ta„hil), (5.7) 

where z^, are the poles of the retarded propagator, and the sum over all possibloe poles is 
taken. As was shown in [16], when a self-energy and hence poles include imaginary part, 
the following substitution of hyperbolic functions with F functions should be made [24] 



'2 cosh(7r2;) ' 

\niz)\' = ■ I, y (5.8) 
z smh(7rzj 

Using the relation between the effective potential and the pressure, p = ^Fe//, we get an 
equation of state 




p= / T^r^ > ; I -JTin (^|r(^ + ^r) + ^ gir;!),. j ■ f^-') 

where summation over complex poles is performed. As in our previous calculations we 
take only the contribution of the nearest to u; = pole eqs. (3. 20), (3. 21), (3. 22), and the 
self-energy T,(z) ~ z'^'^. Near the Fermi surface, the one- loop contribution dominates over 
the self-energy term for the Fermi liquids u > ^, while the sel-energy becomes important 
for the non- Fermi liquids ly < ^. 
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Having calculated the pressure, we can obtain other thermodynamic quantities, e.g. 
the entropy, the specific heat, and the particle number density, respectively, 

with /i = kp- 

We find the temperature dependence for the specific heat. The first term in eq.(5.9) 
contributes to the specific heat 



1 f (Pk „ / , iz* 1, 1, 



, ,^Re\ z^T^i—^ + -); ^ z*T^( * + -) ] (5.11) 

r2 J (27r)2 V 27rr 2' * ^ 2'kT 2' J ^ ^ 

where \I/'(x) = ^ = '^^^"^ ■ The second term in eq.(5.9) gives to the specific heat the 
following contributions 



^ / (0^^' (~ TE{z.)Fn ~ z.E(z,)F[r]; ~ ^^^F[r]) , (5.12) 



where F\r] denotes a combination of T functions and its first and second derivatives. 
Here momentum integration is performed around the Fermi surface, d'^k — t- kpdk± with 
k± = k — kpi the poles = — iV are given by eqs. (3. 20), (3. 21), (3. 22) for the three cases 
of interest, and ^(z) ~ z^'^ . 

For the Fermi liquid > ^, 2;_l ~ /c_l (the real part is dominant). The first term gives 
Yi f dk_\_z1 —7- T and the same behavior from the other combination. In the second term 
we have S ~ k"^^ . Therefore the second term gives ^2 J dk±'E{z^,)z^, — )■ T^'" and the same 
behavior for the other two combinations. In eq.(5.9) for the pressure, the one-loop term 
dominates over the self-energy for u > ^. Therefore for Fermi liquid at low temperatures 
we have 

c~r. (5.13) 

This result reproduces correctly the linear temperature behavior of the heat capacity known 
for the Fermi liquids. 

j_ 

For the non- Fermi liquid < 5, z± ~ kf^ (for both real and imaginary parts). The 

i 1 ^1 

first term gives 1^ f dk±k^ — )• T^"^ and / dk±k^T — )• Tz^. The second term gives 

f dk±T,{z^,)T — )■ r^'^ and subleading behavior for the other two combinations. For v < 
the self-energy dominates over the one-loop in the pressure. Therefore for non-Fermi liquid 
at low temperatures we have 

C'-T^'". (5.14) 

This result for the heat capacity refiects the scaling behavior of the self-energy. 

For I' = ^ all obtained above terms are ~ T. Therefore for the marginal liquids we 
have c ^ T. 
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We repeat derivation of equation of state using the spectral function eq.(5.4). Density 
of states can be written through a spectral function as follows 



/ (02^(^^™'^^) ^ -4^/ ldzA{z,k)fiz), (5.15) 



n 



where f{z) = tanh(2|^). One can use also the Fermi distribution function f{z) = n{z) with 
a prefactor (—2^), which gives the same result for observables. The pressure is given by 

p = dfi'n, (5.16) 
^—00 

where in our case fi = kp- For simplicity we again take hivp — t- 1. We expand the spectral 
function eq.(5.4) with respect to the imaginary part of the self-energy, which we consider 
to be small in this calculation [28] 

A{z,k)^2Tr6{z-z^)-j:2{z,kF)V'- ^ 



z 



where the pole of the propagator is a solution of the equation z — vpk± — T,i{z, kp) = 0, 
which does not contain imaginary part of the self energy S2. Substituting this representa- 
tion in the equation for the pressure, we have 

1 d'^k 1'°° / 1 \ 

" = -4^ L " -'-^^ ^^''>'"— ) <^-'«) 

The frequency integral in the first term gives familiar expression for the number density 

where usually / is a Fermi distribution function, and the dispersion relation is given by 
(in standard notations z^, — t- e^). Here we have f{x) = tanh(|). Therefore integrating 
over kp gives / dk'p tanh^ — )• ln(2cosh^) where to the leading order ~ (/c — kp). 
In the second term we interchange the order of integration in z and kp. Therefore, 
f ^ dk'pV' , / J \ ^ — )• — — nrh — -, to the leading order there is no kp dependence in 
^2(2) ~ z"^^. The second integral is dzTi2{z, kp)f{z)j^ — )■ S2(z*)/(z*). Combin- 

ing all the terms together we have 

f = /(0E(5^>"(2-i'|f) + i^=(^.)'-i'|f). 

where is the pole of the fermion propagator without the imaginary part S2, and there 
is a summation over poles. If we take to be the pole of the full propagator, becomes 
imaginary and generalization of hyperbolic functions to the F functions is necessary eq.(5.8). 
We arrive then at eq.(5.9) for the pressure of the system. 
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5.2 DC conductivity 

We calculate the DC conductivity in the boundary theory using the "dressed" by gravity 
retarded/ advanced fermion propagators eq.(5.3). To make the calculations complete, we 
need the "dressed" vertex, to satisfy Ward identities. As was argued in [22], the boundary 
vertex which is obtained from the bulk one can be approximated by a constant in the 
low temperature limit. Also, according to [26], the vertex has only singularities of the 
product of the Green functions. Therefore, dressing the vertex will not change temperature 
dependence of the DC conductivity [26] . 
Using linear response theory, we have 

a = —^lmUAA{uJ,k = 0)\^=o, (5.21) 
which is a Kubo formula for conductivity. Here the polarization operator is given by 

J2^'^ (^(i^m + iVn, k)KA{iuJm + i^n, ii^m, k)G{iuJrn., k)A^^\k), 

(5.22) 

where the fermion frequency is uim = (2?Ti + 1)ttT, and boson frequency is Vn = 2mrT, 
and in the low temperature limit A^(ia;m + ii^n,i^m,k) = A^^\k). Usually the most 
difficult step is to take the Matsubara sum. Here we do it in two ways. First, analytically 
continuing in the complex plane iujm — ^ z and replacing the Matsubara sum by the contour 
integral with the Fermi distribution function n{x) = whose poles are at the Matsubara 
frequencies along the imaginary axis. Second, using the spectral representation. In both 
cases we follow [26] , where transport coefficients are calculated with propagators including 
imaginary parts. 

In the first way, we have for the fermion Matsubara sum 

H(iun, k) = T'S^ Giiujm + ii^n, k)G(iu!m, k) —?■ -rr^ / dzG(z + iun, k)G(z, k)n(z), 

(5.23) 

where the contour along the imaginary z-axis can be deformed to contour C which goes go 
along the two brunch cuts, ImZ = and Imz = — f„, and the large arcs F with vanishing 
contribution [26]. The fermion propagator has a branch cut along Imz = [27], [26]. 
Therefore we can rewrite 

Hii'^n) = -7^ r dCn(C)G(iz.„ + C)(Gjf(C) -Ga(C)) 



27ri 
1 
2Tri 



oo 



dCn{C)G{-iu^ + OiGRiC) - Ga(C)), (5.24) 



where the difference of the retarded and advanced functions in the first bracket is due to 
the discontinuity along Imz = and in the second bracket due to the discontinuity along 
Imz = —Vn- This contribution corresponds to the second term in eq.(5.23), and there are 
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no pole contributions [26]. We use the usual prescription for the retarded and advanced 
Green functions, Gr = G(a; + iO"'') and Ga = G{oj — iQ'^). We suppress momentum indices. 
Taking ivn — )• cj + iO"*", we have 

= dQn{C)GR{u + 0{GR{0-GA{Q)) 

- ^ r dCniC + uj)GAiu; + C)(Gr(C + GAiC + ^)), (5.25) 

where we changed the integration variable in the second integral (—uj — t- In the limit oj — )• 
0, the dominant contribution comes from the pair GrGa, and it is inversely proportional 
to the distance between the poles given by the imaginary part T,2- Combinations GrGr 
and GaGa with the poles on one side of real axis make a much smaller contribution due 
to cancellation between the residues at the poles. Therefore, as ~ 0, we have 

H{io, k) ^ --^ r dCiniC + co)- niC))GR{C + co)Ga{C), (5-26) 
2Tn 



lmUAA{u^,0) = ^ I -^Afik) r ^(n(C + a;)-n(C))Gfi(C + ^,fc) x 



and 



2nJ (2vr)2-^ ^'^V^oo 2vr^ 
Aa(C + 0J + iO+, C - iO^,k)GA{C, k). (5.27) 

In small T limit the vertex is a constant. We integrate around the Fermi surface, therefore 
momentum integral is J — ^ ^(2!!)^ with k± = k — kp. We exchange the order of 
integration and do first the momentum integral [14], [22]. For w ~ 0, we have 

dk± 1 



2^ (4 - fc± + S(C, kp) + i0+)(4 -k^ + S*(C, kp) - iO+) 
i 1 



s(c,M-5:*(C,M 2imS(c,M 

Writing n'(C) = — /3n(C)(l — n{C)), we have for (x) ~ 



(5.28) 



where we did not include constants. Note that we get the same result for conductivity 
when we use tanh | in the contour integral eq.(5.23) since n'{x) = — 2tanh'(|). 
For the Landau Fermi hquid T,{uj) ~ uj^ at small T [29], [22]. We get 



a ~ T-^ (5.30) 

that means we recover the standard result for the resistivity of the Fermi liquid, p ~ T^. 
In our case, 5](tj) ~ uP''^^f ^ 



a 



T-^^^F^ (5.31) 
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where Vkp^s the IR conformal dimension. This result agrees with the DC conductivity 
obtained in [22]. For the marginal liquid, u^p = ^, we recover the resistivity p ^ T, 
which is known for strange metals. It is interesting that the scaling behavior of the DC 
conductivity is the same as the single particle scattering rate. In the gravity calculations 
it is explained by the fact that the dissipative part of the current-current correlator is 
controlled by the rate of the bulk fermion falling in the horizon, given by the single particle 
scattering rate. 

To check our calculation, we get the DC conductivity using the spectral representation 

G{iUm,k)= I —- , (5.32) 

J 2tT ko- lUJm 

where the spectral function A{ko,k) is given in eq.(5.4). For the product of the Green 
functions we use the following formula 

^ ioJm. - Wl iuJm + il^n - ^2 i'^n + UJi — U}2 
m 

Taking ivn — )• w + iO^, the polarization operator is given by 

^ / N f d?k dui dujo n(uji) — n(uj2) , (0)2 , , , n w , n 

n.x{-.0) = J p^^^ ^ a!?' A(.„k,)A(^M. (5.34) 

Performing integration over 002, we have 

lmUAA{uJ,0) = J ^;^^{n{u;i)-n{i02))I^TM^uk±)A{oJi + u;,k±). (5.35) 

In the limit a; ~ 0, the momentum integration is 

with S2 = ImS. Therefore the DC conductivity given by eq.(5.21) is 

^ V 27r lmJ:{uJi,kF) ^ ' 

which is the same as eq.(5.29) obtained by the contour integration. 

6. Discussion 

In this article we studied the particle-hole pairing in the context of the magnetic catalysis. 
The Reissner-Nordstrom charged black hole can carry Fermi surfaces [14], in a sense that 
one has a Fermi liquid in the bulk: there are Fermi surfaces and free fermions fill all the 
levels up to the Fermi surface. It is natural to expect pairing between fermions (particle and 
hole in this case) when an attractive interaction is introduced. At the CFT boundary, we 
indeed get pairing instability for the Fermi liquids. However, quite surprisingly, there is no 
pairing realized for the non-Fermi liquids. We show that in variational calculations where 
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both terms in the r.h.s. of the gap equation vanish near the Fermi surface. In the bulk 
Ginsburg-Landau calculations, we obtain the critical temperature which vanish exactly at 
the IR CFT conformal dimension V]^p = ^, which is the border between Fermi and non- 
Fermi liquids. Tc stays zero for < ^. The same conclusion has been reached for the 
superconducting instability in [14]. Probably the reason is that there are no well defined 
quasiparticles in case of non-Fermi liquids, in particular the residue for the fermion pole 
vanishes around the Fermi surface [18]. Note, that we have well defined particle degrees of 
freedom in the bulk, while the situation is different when it is projected to the boundary. 
As suggested in [14], the momentum/frequency dependent four- Fermi interaction may help 
to realize pairing for the non-Fermi liquids. 

Our calculations are different from the corresponding calculations in the field theory 
in several aspects. We have imagined the existance of an "experimental knob" which can 
be used to adjust the UV scaling dimension A of an operator, and therefore an ability to 
get different CFT^s, e.g., describing Fermi and non- Fermi liquids. As discussed in [18], the 
most useful knob will likely depend on the UV geometry into which this AdS2 is embedded. 
For example, in our case an external magnetic field will allow one to tune the IR scaling 
dimension u. We obtain a radial profile for the order parameter A(r) ~ '(/;'^(r)V^'i/''^(r) 
where ip'^^r) are the zero modes with oj = and k = kp. For the Fermi liquids, Ukp > 5, 
A(r) is supported near the CFT boundary, and for the non-Fermi liquids, Ukp < |, A(r) 
is supported near the horizon where the Fermi velocity is considerably smaller than the 
speed of light vp <^ c. The radial profile of A(r) with correct fall off is important to insure 
convergence of the radial integrals. 

In the presence of the particle- hole condensate < ■^p^T■^T^^p >, the Fermi momenta 
corresponding to the bulk fermions with spin up and down are shifted in the opposite 
directions. We did not introduce spins in the boundary theory. Therefore association with 
the bulk spins is understood as upper and lower components of the 4-component fermion 
field. For a fixed fermion charge we obtain multiple Fermi surfaces. This is the consequence 
of the AdS geometry effectively behaving as a box potential. It is interesting to understand 
the physical picture behind the multiple Fermi surfaces, and if pairing between different 
Fermi surfaces brings new physics. 

In this paper the particle-hole pairing was considered on a very general ground. There- 
fore its evidence can be relevant for the chiral spirals [3] in the chiral magnetic effect [30] 
and the spin density waves, and can serve as a guide to construct an antiferromagnetic 
and Mott insulating states of the cuprate superconductors. In particular, it can be use- 
ful to describe the coexisting AFM and SC order parameters in the iron pnictides [31]. 
Recent work on non-abelian holographic superconductors [9] makes applications to color 
superconductivity possible. 

It will be instructive to obtain, in analogy to the superconducting instability, to obtain 
the particle-hole condensation from the classical Einstein gravity calculations with the 
neutral scalar field. 

We considered thermodynamic and transport properties of the boundary field theory. 
We did not specify the boundary theory. Using the "dressed" by the gravity fermion 
propagator Gr we obtained the equation of state, which includes a sum over all the poles 
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of Gr. We use prescription of [16], [17] to treat imaginary poles. Imaginary self-energy is 
a consequence of inelastic scattering by the black hole, and will be present in any gravity 
calculations. In particular the imaginary part provides finite transport coefficients. The 
scaling behavior for the heat capacity is c ~ T for the Fermi liquid and c ~ T^'^, while it 
is the same in both cases for the DC conductivity a ~ T"^'". It reflects the fact that the 
specific heat depends on the dispersion relation, in particular on the real part of w, while 
the DC conductivity is sensitive only to the scaling of the self-energy. 

The presented approach has an advantage of unifying description of the Fermi and non- 
Frmi liquids within one framework. It uses the language of poles and branch cuts instead of 
quasiparticles, that may be more adequate description for some strongly correlated systems. 
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A. Dirac equation in the AdS^ 

Here we discuss the Dirac equation in the presence of the magnetic field the AdS^, and show 
how Landau levels appear for the {x, y) part, dimensions of the boundar field theory. The 
(x, y) part of the Dirac equation decouples from the radial part and can be solved exactly 
due to the translational invariance in perpendicular directions. Then we find the confor- 
mal dimension of the spinor operator in the IR CFT. Depending on the IR conformal 
dimension the boundary theory describes Fermi, marginal and non- Fermi liquids. 

A.l Dirac equation with magnetic field in a charged black hole geometry AdS4 

Here we solve analytically the part of the Dirac which depends on magnetic field and space- 
time coordinates of the boundary theory. The free spinor action in the geometry eq.(2.2) 
and in the presence of magnetic field eq.(2.3) is given by 

So = iJ d^xV^^P {T^'Vm - m) V^, (A.l) 

where ■0 = i/j^T^, and 

'Dm = Om + ^uJabMT"'' - iqAM, (A.2) 

with ujabM the spin connection, and F'*'' = |[r'*, F^]; here M and a, b denote the bulk space- 
time and tangent space indices respectively, and fi, v denote indices along the boundary 
directions, i.e. M = {r,n); indices with hat refer to tangent space ones, i.e. converting 
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metric eq.(2.2), ds^ = guNdx^ dx^ = Vah^"'^^ ^db ~ '^i^sl"!; 1; 1; 1) metric. 



from bulk to tangent indices = e^F"' with e° = e'^.jdx^^ are the tetrads defined by the 
xic eq.(2.2), ds^ = guNdx^ dx^ = 
Using the translational invariance, 

il}{t,x,y,r)= I dujdke-''^^+''''' 4j{u},k,y,r), 



(A.3) 



with k = kx, the Dirac equation is given by 



1 



1 



1 



1 



T'i-iu + -Ui^^T'^ - iqAtir)) + -—T'dr + —T'^iik + -oo^f^T^' - iqA,{y)) 



+ -^Ty{dy + l:UJyfyryn - ij{uj, k, y, r) = 0, 

V 9ii ^ J 



'9v. 



(A.4) 



where gu = g^x = 9yy, and At{r) = fi{l - ro/r), A^iy) = -Uy. From the torsion-free 
condition, A e'' = -de"", we find the spin connection [32] for the metric 2.2, 



tf 



^= — dt, uj:~ = dx , 



where i = x,y. Note that 



and 



e^T°'uj' r 



(A.5) 



(A.6) 



be 



4 « 



be A'/-" 



r + 



1 1 dr{y/- gtt)^f , 2 1 dr^i 

4 V-S'it a/^^ 



1 



In 



1/4 



where g is the determinant of the metric. Therefore, we can rescale the spinor field 

-1/4 



9_ 

9rr 



and remove the spin connection completely. The new action is given by 



5n 



(A.7) 



(A.9) 



where the covariant derivative does not contain spin connection, V'^ = Om — iqAm- 
In new field variables, the Dirac equation is given by 



r'dr 



/gr 

+ Vydy) <^iu},k,y,r)=0. 
We separate y and r dependences, 

P(r) = -^/^T^dr 



. =r* i(a; + - ^)) - + i{k + qUy) 

v-gtt r 



(A.IO) 



r + 1 ) 

r 



grr \/—gtt 

Q{y) = V"^ i{k + qny)+Tydy, 



/gum, 



(A.U) 
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and the Dirac equation is 

{P{r) + Q{y))'^ = 0. (A.12) 

Though, [P{r),Q(y)] ^ 0, one can find a transformation matrix U such that [UP{r), UQ{y)] ■ 
0, and then look for common eigenvectors of UP{r) and UQ{y) since they are commuting 
hermitian operators, i.e. the Dirac equation reads 



C/P(r)$„ = -UQ{y)^n = A„,$„, 



(A.13) 



where n will label the Landau levels. In the main text we use I for the Landau index, 
in order not to confuse with the Matsubara frequency index n. Transformation matrix U 
should satisfy conditions 



{u, r } = 0, {u, r*} = 0, [u, r^i = o, [u, ry] = o, 

which do not fix U completely. It is convenient to use the following basis [18], 

^ .3 ^ 



(A.14) 



—a 
- 

■ <J^ 




-a 
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icj2 
-icj2 



(A.15) 



Note, that the usual relation holds 

p5 p0pl"p2p3 

with —7- 1, 1 —7- X, 2 —7- y, 3 —7- r. In the representation of eq.(A.15), we can choose 

" \ 



(A.16) 



U 







-la 



(A.17) 



Writing $ = (Fi, Fs)^, and using eq.(A.31), we get the Dirac equation written in a compact 
form, eq.(A.13), 



'9ii 1 o , , . 2 

<T Or + y/guKJ m 



V-9tt 



a^{uj + qfi{l - ro/r)) - An 



1 



-{k + qHy)+\n 
-id,, 



id,. 



+ qHy) + A 




0, 




where m. X®Y , X acts inside Fi or F2 and Y acts between Fi and F2. In eq.(??), the 1 in 
the first equation shows that there is no mixing of Fi and F2 by the operator UP{r) and 
the 1 in the second equation means that there is no mixing inside Fi or F2 by the operator 
UQ{y). Therefore solution can be represented as 




(A.20) 
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we do not write uj and k dependences. Dirac equations for each component are 



/gv 



9rr 
\f9r 



dr^^rn] /«(r) + 



I9i 



\/-9tt 



(oj + qfi{l - ro/r)) + A J )(r) = 0, 



'^mj fPir) + (^^(^ + - ro/r)) + \^ A'\r) = 0, 



(A.21) 



-idyg^nHy) + {{k + qHy) + A„) ^i'^ = 0, 
-idyg^^^ + {-{k + qUy) + Xn) gl^^ = 0. (A.22) 

In equations A.22 for the y dependence, we rescale y = \/\QTT\ {y + k/qH) and An = 
^/\qT-L\ Xn, and get 



-idyg^n^ + {y + An)*?!') = 0, 
-idyg^^) + {-y + Xn)g^^^ = 0. 



(A.23) 



The second order ODE 



1 



y + x 
1 

-y + x 



^dyg('^ + (A' - 2/2)^(1) = 0, 



dyg^'^ + {X' - yV^ = 0, 



(A.24) 



are 



are solved by substitution g^^^ = e~^'^^'^g^^^ and gf*-^) = ■+.iQ~^'^ The eigenfunctions 
Hermite polynomials. We get the same eigenvalues, but slightly different eigenfunctions 
for different signs of qH.. Putting all together, for qK > 0, we have 



A_i = : 

Xt = ±\/2(n+l) : 



(A.25) 



V y'2{n+l) 



Hn+i{y) 



Hn+i{y) 



and for qH < 0, we have 



A„i = : 

/\)(y)=e-^^/^ g^^\{y) = ^.-y'l' 
Xi = ±^2{n+l) : 

1 



(A.26) 



g-. 



V2(n+1) 



Hn+i{y) 



g^^Hy) = -i^-''" {^Hn{y)±-^ 



+ 1) 



Hn+iiy) 
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The case qH < can be obtained from the case qV. > by replacing g^^^lqV. < 0] = 
-ig^'^^'iqn > 0)] and g^^'>[qn < 0] = -ig^^^qU > 0]. Using the eigenvalues, eqs.(A.26,A.27), 
in the equation for the radial part, eq.(A.21), we get 



'Wr J 

+ f-^iL(u; + q^l{l - roM) ± V2\q'H\in + l)) = 0, 

V V-9tt ) 



+ + 5^(1 _ ro/r)) ± ^2\qn\{n + 1) ) /«±(r) = 0. (A.27) 

\\J-9tt 

These equations can be obtained by replacing 



k ±V2|g^|(n + 1), (A.28) 

in the Dirac equation at zero magnetic field [16]. Equation (A.28) also gives a prescription 
how to treat the limit of zero magnetic field, i.e. the — )• limit is taken keeping 
2\qH.\{n + 1) = A;^ fixed as — )• 0. In a compact form the Dirac equation, eq.(A.18), in a 
magnetic field reads 

( ]=&^dr - m H — (j^(a; + g/i(l - ro/r)) 

V \J-9tt 

^ J__ia^ ^2\qn\{n + 1)^ ^ 1 ^ ^ = 0, (A.29) 

which coincides with eq. (A14) in [18] with the replacement eq.(A.28). The starting Dirac 
equation in a magnetic field, eq.(A.lO), is given by 



-r'dr - ^=r* i{uj + qi^i{l - -)) - m 



^ V-9tt r 

T -^U^W2\qn\{n + 1)] <^{r) = 0, (A.30) 
V9ii J 

we do not write uj dependence, $ = {Fi,F2)'^, n = —1,0,1,...; and where is the 
inverse matrix to the transformation matrix eq.(A.31) 

We use equation (A.30) in the main text. 

A. 2 Dirac equation. Conformal dimension in the low frequency limit 

As outlined in [16, 17], the fermion determinant in the black hole background is given by a 
sum over the quasinormal frequencies = iun, that are obtained as complex frequencies 
which give zero eigenvalues, A = 0, of the Dirac equation. The eigenvalues are given by 

(t^^Vm -m- AiT^r^ ± qnr^^ V = AV'- (A.32) 
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In Appendix A.l, we simplified the Dirac equation obtained from tlie free fermion action 
eq.(3.2). We rescaled tlie original fermion field, eq.(A.8), as was done in [18], and removed 
the spin connection, with the result eq.(A.9) and eq.(A.lO). In new field variables, the 
Dirac equation is 

r^^v'j^.^ -m- Air^r^ ± qnr^) $ = o, (A.33) 



where T>m = Qm — ^qAm- It is written in the geometry eq.(2.2) and in the Landau gauge 
as 



-LrS, - —^iT\u; + g/i(l - !^)) - m + -^iT^ik + qTiy) + -^r^dy 
iQ^ V~9tt r WgTi Jg^i 



with k = kx, ky = 0. At the boundary, r — t- oo, the frequency oj is measured from the 
effective chemical potential qfi. The dependence on charge q of tp-field enters only through 
the combination qfi. Also as given in Appendix A, in the Landau gauge eq.(2.9) the radial 
and y dependences decouple, eq.(A.18) and eq.(A.19), and y dependence reduces to the 
harmonic oscillator. We solve the y dependent part of the Dirac equation, eqs.(A.26) and 
(A. 27), and use the result in the radial part of the Dirac equation. The procedure of 
obtaining the Dirac equation at nonzero magnetic field amounts to replacing 

k ±y'2\qn\l (A.35) 

in the Dirac equation at zero magnetic field. As compared to eq.(A.28), eq.(A.35) takes 
into account the Zeeman splitting term. The result for non-interacting fermions is given 
by eq.(A.30), which agrees with [16, 17]. The Dirac equation is written as 

^rdr - —^iT\uj + g//(l - ^)) - m ^ -^U^^^/2\^l - AiV^T^) $ = 0, 
V-9tt r y/gTi J 

(A.36) 

where ^ = 0, 1, . . .. The F matrices are defined in eq.(A.15) as in [18] and is given by 
eq.(A.31). 

We split the 4-component spinors into two 2-component spinors (we do not write zero 
entries) F = {Fi, F2)^ where the index a = 1, 2 is the Dirac index of the boundary theory, 
using projectors 

n„ = ^(i-(-i)°rW), a = 1,2, ni+n2 = i, (A.37) 

which commute with the Dirac operator of eq.(A.42), and Fa = IIq.^, a = 1,2, decouple 
from each other. Gamma matrices were chosen in such a way that this decoupling is 
possible. The Dirac equation for two components is given by 



^dr-m+ -^=(T^{u: + q^i{l - ro/r)) T -^ia"^ ^ 2\qn\l + (-l)"sA ) F^ = 0, 



1 

a"Or — m , ^„ , ^^y^ -u/'/y I 

V-9tt ^^9. 



(A.38) 
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where s = isgn{q'H). To obtain the retarded Green function for fermionic operator O in 
the boundary theory, we need to find a solution ^> which satisfies the ingoing boundary 
conditions at the horizon, and to expand it near the boundary at r ^ oo as 



F2 « aar"^^^ |^ J ^ + ^sr"™^^ |^ J ) , (A.39) 

where is the source and 6^ is the v.e.v. in the Green function terminology, and the two 
masses are mi = m + A, m-2 = |m — A|. The conformal dimension of O is given in 
terms of the mass of i/j field by eq.(3.1), A^ = i + nT-iR- The two spinors (1, 0) and (0, 1) 
are eigenspinors of with opposite eigenvalues, implying that Cq, and 6q are canonically 
conjugate (in a radial Hamiltonian slicing) [23]. Depending on the value of the exponent, 
we impose different boundary condition for and/or ba- The retarded Green function is 
given by [33, 23, 34] 

GRiu, k) = ( ^) ° ] , with Go^iu;, k) = ^, a = 1, 2, (A.40) 

\^ G2{uJ,k) J Qa 

with k y/2\qn\l. 

For a general Reissner-Nordstrom black hole background there is no mixing between 
Fi and F2 (interaction term Sint is diagonal), and in the absence of A-term two fields Fi 
and F2 have coincident Fermi surfaces (at a; + g/i = at the boundary r — t- 00). With the 
A term, there is a relative shift of the Fermi surface in the spectrum of Fi and F2. In the 
following, we consider only one component, e.g., Fi. 

We could have solved eq.(A.38) with the ingoing boundary conditions at the horizon 
numerically. Instead, we use analytical results of [18], and for that we need to extract 
information from the low-energy limit of the theory. In section 4, we follow [14] and solve 
the Dirac equation at zero frequency, and then use the matching procedure to extract the 
retarded Green function of the boundary theory. For now we use results from [18, 20]. 
We will be interested in the poles of the retarded Green function. As was shown in [18], 
the nonanalytic part of the retarded Green function in the extremal Reissner-Nordstrom 
black hole background eq.(2.2) comes from the IR region, which has a simpler geometry 
than AdS4. At small frequencies, — t- 0, the metric reduces to AdS^ — t- AdS2 x R2 at the 
horizon 

ds^ = ^i-dT' + de) + '^dx^ Ar = ^^- A, = -ny, (A.41) 
which is called the IR region [18]. In the IR region, the Dirac equation (A. 38) becomes 



=-'^C - m + ^cx^(- + ^4 J) - - «a) F = 0, 

5CC v—9tt v12'^**<, j 

(A.42) 
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0. (A.43) 



where we omit index by the spinor field. We can write exphcitly 

At the AdS2 boundary, C ~^ 0) the Dirac equation to the leading order is 

Ca.^ = U=\ ' ^ ''I (A.44) 



Diagonalizing matrix [/, we obtain the conformal dimension for the operator O in the 
IR CFT dual to 

1 

<>V' = 2 + ^' 



(m + .A)2 + "^215^1^ Ui - (A.45) 

where the AdS^ curvature radius is R2 = R/VQ, the Landau level index I = 0,1,.. . and 
s = ±sign{q'H). The conformal dimension of the fermionic operator in the IR CFT in the 
chiral limit m = and at the Fermi surface is given by 



'(M^ + A2)i?i-^^, (A.46) 



with kp = qfi, /i = \/3(7f;^^. We use the IR conformal dimension of the fermionic 
operator in the main text. 



B. Dirac equation in the AdS2 

Here we discuss the Dirac equation in the AdS2- We obtain the conformal dimension of 
the spinor operator in the CFTi, which gives (with trivial modifications) the conformal 
dimension of the operator in the original CFT2-\.i in the IR region. It controlls the behavior 
of the theory, e.g., Fermi liquid or non-Fermi liquid regimes. Then, we derive analytically 
the two point Green function, which up to a nonimportant constant defines the self-energy 
in the fermion propagator of the CFT2+1. 

B.l Dirac equation and conformal dimension 

We consider the following action for a two-component spinor field i/' 

S = iJ d^x^^i) [v'^Da -m + iT\- FA) V', (B.l) 
with Da = da — iqAa and the background metric and gauge field given by 

ds' = §{-dr^ + de), ^r = ^\. (B.2) 
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In eq.(B.l), we included a time-reversal violating A term which in our application will be 
related to kx momentum in R^, and hence the eigenvalue incorporating effect of the 
magnetic field; we also included A term to mimic the gap in R^. We choose the following 
Gamma matrices 

^^ = -cJ^ r^ = ia^, r = -a^, f = i, (b.3) 

where the hat indices denote those in tangent frame and a* are standard sigma matrices. 
Writing 'i/'(r, (") = {—g/gc_(^)~^^^ J du}e~^^^{u}, C) the equations of motion are 

—a^dt - m„ + , ^ a^juj + qAr) - ia^x] ^> = 0, (B.4) 
/5CC V-frr / 

where we introduced m** = m + A. In a matrix form we have 

To find the conformal dimension of the operator O dual to ^, we solve the Dirac equation 
near the boundary. Near the AdS2 boundary (" — )• 0, equation (B.5) to the leading order is 




0. (B.5) 



As C — 0, can be written to the leading order in the form 

^ = Av+C"" + Bv-C, (B.7) 
where v± are real eigenvectors of U with eigenvalues itz/ respectively 



v = ^{{m + A)2 + A2) _ (B.8) 

with R2 = R/VQ. Here u determines the scaling dimension of the fermionic operator 
in the CFTi dual to the spinor field in the AdS2- With small modifications this gives 
the scaling dimension of the CFT-^ in the IR region, since as w — t- the metric reduces 
AdSii AdS2 X R2 near the AdS^ horizon. The eigenvectors of U can be chosen from 
eq.(B.6) as 

/ (m + ^)R2 , ^ 



n2rl, / 

Imposing the ingoing boundary condition for $ at the AdS2 horizon, the retarded Green 
function for the boundary operator in the IR CFTi dual to <1> can be written as 

GW{^) = f a;2^ (B.IO) 
meaning that the operator conformal dimension is given by 

<5^ = ^ + i/, (B.ll) 

where u is the scaling exponent with respect to w, given by eq.(B.8). In eq.(B.lO), we used 
the fact that in the solution oj and C, scale with the same power [18]. 
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B.2 Two-point functions for charged fermions in AdS2 

We calculate G^^ analytically. The equation of motion (B.4) can be written as 

= ia^ L + ^f^-^ $ - ^ (a^m.. + a'\) (B.12) 



where m** = m + A. The equation of motion for two components $ = (y, z)^ is given by 

2 

(C^c + m„i?2)y - K + - Ai?2)2 = 0, 

2 

(C^c - m„i?2)z + (wC + + Ai?2)2/ = 0, (B.13) 



which contain C, dependence in the mixing term proportional to io"^, that makes it difficult to 
solve for one component. We therefore transform to another basis to make all C, dependent 
terms diagonal [?]. We make the following basis rotation 

Transforming the sigma matrices, M(T*M~^, we have ia"^ — )■ ia^ , — )• — cr^, cr^ — )• cr^. The 
equation of motion becomes 

ac4> = ia'(u: + ^^4-^1 + ^ (a^m,, - a^\) ^, (B.15) 
and for two components ^> = (y, z)"'", 

2 

(C5^ - i{ujC + ^E^)) y + (im,. + A)i?2 5 = 0, 
V 12rf 



.2 



{Cd. + i(wC + ^^^)) z - - A)i?2 y = 0. (B.16) 
V 12rf 



Expressing y from the second equation 



y = (B.17) 

we get the equation for z. The equations for both components are given by 



Rewriting these equations, we get 

^2 



^ ^ ^ ^S'Fr^ x\ - / 2 , ^2^D2- 



12r2 



^ = (mt + A')i2^5. (B.18) 



eOlij + C^cy + (-^Cc^ + (a;C + ^f^f - (ml + X')RI) y = 0, 
C^a^g + ^Q^~, + + (^^ + ^f=4-)' - + = 0. (B.19) 
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MATHEMATICA gives the following solutions for equations (B.19) 



m 



-iuiC, /-V 



l + 2i/, liujQ +C2L{-1 



z{C) = e-'^'^C ( C3 U{u + 1 + 2u, 2iojC) + C4 L{-v - i^^£^,2u, 2iojC)] , (B.20) 

* J 
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where u = \ {"ml^ + A^)i?2 ~ '^127^4'"* 1 U{a,b,z) is the tricomi confluent hypergeometric 

function (of the second kind) and L(z^, A, z) = L'^{z) is the generalized Laguerre function 
(for V = n the associated Laguerre polynomial). We substitute solutions (B.20) into the 
system of first order ODE, eq.(B.16), and consider this system at the AdS2 boundary, 
C — 7" 0, where it is considerably simplified. As a result we get the relations between the 
constants 

ci 



— = (im** + A)i?2 

C3 



C2 
C4 



{im^^ + A)i?2 



I2rl 



{—im^^: + A)i?2 ' 



(B.21) 



to simplify we use r(z)r(l — z) = ir/ s'm{Trz). Both equations in the system give the same 
relations eq.(B.21). As a consistency check, we found the same relations considering the 
system at the AdS2 horizon, ^ — >• 00. In order to insure the ingoing wave ~ e*"^^ at the 
horizon — )• 00 for each of the solutions, we have one more relation between the constants 



ci 

C2 



12r2 



TT 



(B.22) 



Relations (B.21, B.22) fix all the constants up to overall normalization constant. Using the 
relations (B.21, B.22), the solution of the system (B.16) at C = becomes 



{2iuj) 



{2iuj) 



-2v 



-2v 



r(2i/)r(i 



12rL 



r(2z^)r(i 



12rL 



T{-2v)T{l + v 



12r2 



-V + i 



12rL 



12r2 



r{-2iy)T{l + I' 



12ri 



At C 



(im** + A)i?2 (im** + A)i?2 

0, solution to the original system of equations (B.13) is given by 



-C, (B.23) 
I 

-C- (B.24) 




{2iu}) 



-2u 



T{2v)Y{l 



12rL 



2{v + i 



12r2^ • 



13 pr 



l + i 



\ 



{ims,s,+\)R2 



\ {im.„+X)R2 + * / 



+ 



r(-2.)r(i + .-z^ 



2(-u + 7 ) 



1 + ^71 



qspr 



\ 



(imtt+X)R2 
2 



V (I 



-I/+1 



13 pr 



+ 1 



(B.25) 
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The system of equations (B.13) near the AdS2 boundary has a solution which can be 
written in a general form as, eq.(B.7), 



y 



then the Green function is G^j^{u}) = B/A. We put solution eq.(B.25) into this form, 

v+i 

^2 



z 



(2zt^) r(2z.)r(i-z/-z-=— )- — — -\ ^g^rl 



12ri, 

+ r(-2.)r(i + . - z ^ ' + ^^^^^ ' ""''^ " ' 




and extract the IR Green function to be 



V{-2v)V{l + u-i^^SFTL 



((m„ - iX)R2 + - + v){m^^R2 + v) 

^i^^^ii {2uf\ (B.28) 

{{m..-^X)R2 + ^^^+v){^^^-u){mR,-v)' 



Simplifying the following ratio 



we get the retarded IR Green function given by 



(B.29) 



V{-2v)T{l + u- i^ff^) ((m** + i\)R2 - %^f=r^ - v) 



G^iXuj) = e-'^" {2ojf\ (B.30) 

n2u)T{l-.-^^) {{m.. + ^X)R2-^^^+v)' 

with m=n, = m + A. In eq.(B.30), we did not include a ratio ™**^^^^ , since there is 
an ambiguity in definition of the Green function up to a real function of A (or k with no 
magnetic field) and g, m. If the matching is done using our basis then this difference should 
not matter. This expression for the IR Green function agrees with the one obtained in [18]. 

C. One-loop calculations in a (2 + 1) dimensional field theory 

We calculate here the free fermion energy and the gap equation. One-loop fermion effective 
action in the chiral limit, m = 0, is given by 

5^/^ = -ilndet(i^ - A) = -^lndet(f ^ + A^), (C.l) 
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where ip = {idt + n)^^ — vpK^, and K = iV + qA. For simplicity, we added to the free 
part the interaction Ginti'4'tp)i'4''4') — >• {^{i'lp) + h.c.) — A'^/4Gint, where the strength of 
intercation in (2 + l)-d is Gint ~ ]g^- Here, the order parameter is A = 2Gint < V'V' >■ Iii 
the Landau gauge A = (— ?^y,0), and after the Fourier transform, we have 

-P^ = ioj + fif -vlK^ -iq-UvW^^ (C.2) 

To calculate the fermion determinant, eq.(C.l), we use IndetG"^ = TrlnG"^. The eigen- 
values of operator are known (2/ + (we also calculated them in Appendix A.l); 
the eigenvalues of operator ivpqU'y^^'^ are ±Vp\q7i\ (in standard representation for 7 ma- 
trices); i.e., the H dependent part is (2/ + l)vp\q'H\ ±Vp\q7i\. One can rescale / — Z — 1 for 
one of the signs and combine two terms with both signs together with the result Vp2\q'H\l. 
After rescaling there will be however different prefactors for two signs from taking matrix 
elements under the trace, Tr (see [1] for details). Since we will consider only the lowest 
Landau level, we can ignore the difference in prefactors, and moreover 

(27r)2 ^ (27r) ' ^^"^^ 

that takes into account the degeneracy of Landau levels, since the Dirac equation eigenvalue 
A and hence the quasiparticle spectrum do not depend on momentum k. Here V2 = L^x Ly 
is the size of the sample. We therefore have 

n 

where the fermionic Matsubara frequencies at temperature T are Un = (2?^ + l)7rT (we 



changed to Matsubara frequency by Wick rotation Un = iw), and Ei = ^ 2vp\q'H\l + A^. 
Sum over the Landau levels I is implied. The Dirac equation eigenvalue is A = ((x>„ + i/i)2 + 
Ef, which gives quasiparticle poles z*(/) = iuJn at A = equal to z*(/) = fJ,±Ei. We rewrite 
the Matsubara sum as a contour integral 

, {ujn + ifJ-f + Ef _ i f dz -{z - fif + Ef z 



n 



In ^ " ' — '- = - — In^ — ^tanh— , (C.5) 

^ r2 2 Jc 27r r2 2r' ^ ' 



due to the fact that the poles of tanh are situated along the imaginary axis at z = i{2n + 
1)ttT. Differentiating both sides with respect to Ei, we take the r.h.s. integral 

V , = - E tanh MM = 1 f tanh ^ + tanh ^) . (C.6) 

^ {ujn + ifJ-Y + Ef 2^ 2T 2\ 2T 2T J ^ ^ 



Integrating back over Ei, we have 

\2 I Tp2 



TEln '"" + g = rEln(2coshMi) =r5: (^f + ln(l + eH-™'-)) 

= +rin(l + e-(^'-^)/^) + +rin(l +e-(^'+'^)/^). 

(C.7) 
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A useful formula following from eq.(C.6), 

^K + v)2 + Ef 2Elfj/,2 2r 2E, cosh a + cosh a ' ' 

Putting all together, an effective action for A is given by 




T\<,n\ 



Eln(2cosh5fl)). 



^e// = ^li7=^-i^2^lM2cosh^) |, (C.9) 

with z4l) = fi± El, El = v^2|g^|/ + A2, sum over the Landau levels / is implied. The 
free fermion energy can be obtained by dividing Seff by the space-time volume, i.e., = 
—Seff/{TV2). Minimizing effective action, 6Seff/SA = 0, we get the gap equation 

^^G^A ^nhf 

El cosh f - + cosh ^ 

with El = ^2|g-H|/ + A2, sum over / is implied. For the lowest Landau level, / = 0, the 
gap equation reads 

A = ^^ . (C.ll) 

vr cosh Y + cosh ^ 

At T = 0, the solution is given by 

1 



vr 



Gint\qn\, (C.12) 



provided A > /i, and where Gint = jgp- At T / 0, from eq.(D.2), there is the second 
solution A = 0, and the phase transition between A 7^ and A = 0. The character of the 
phase transition, first or second order depends on the values of parameters [1]. 

We calculate the critical temperature of the phase transition. We fix the charge density, 
n, and express the chemical potential through n. From the effective action eq.(C.9), the 
charge density and the gap equation for the lowest Landau level are given by 



\qn\ sinh( 

n — 



27r cosh(^) + cosh(^; 



_ Gint\qn\ sinh(^ 

We introduce the filling factor 
then from the expression for the charge density, we have 



A = ' ; . (C.13) 

vr cosh(^) + cosh(i^) ^ ' 



27rn Tic ,^ ^ 



cosh(^^) = ■ (C.15) 

1 ^ 



rj^ cosh(f ) + Jl+vli sinh2(f 
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Therefore the gap equation becomes 



^^ G^nMn (l-r7|,)smh(^) ^^^^^^ 



At T=0, the solution is given by 

A=^^(l-r?«). (C.17) 
vr 

There is no nonzero gap for the filhng factor 77-^ > 1. The condition r/-^ < 1 to have a 
nonzero gap translates for the charge density to be smaller than critical one, n < Uc, with 
ric = Ti{r]-}{ = 1), or for the magnetic field to be larger than the critical one, Ti > Tic- 
For rj^ > 1, i.e. n > ric or % < Tic the symmetry is restrored, A = 0. Around the 
critical temperature, when the gap is vanishing, the gap equation gives the following critical 
temperature 

Te = ^^(l-r?U (C.18) 

where Tc = for 77-^ > 1, i.e. for Ti < Tic- For H > Tic, Tc grows linearly with magnetic 
field, Tc ~ Iq'Hl, in the vicinity of the phase transition. Away from the phase transition 
one should solve the following gap equation for the lowest Landau level numerically 

2rcsinh(^) 

A = =. (C.19) 

cosh(f ) + Jl + T]^ sinh2(^) 

We use this procedure to derive the gap equation and to calculate Tc in the AdS4. 



D. Critical temperature from the ADS4 variational calculations 



We calculte the critical temperature Tc for the case of Ladau Fermi liquid, u^p > \. Let us 

introduce an analog of the charge density in the AdSi^ by differentiating an effective action 

s s 

eq.(3.23) with respect to the Fermi momentum kp, n{r) = ^(^^{j^-) ■ Together with the gap 
equation, = 0, we have 



n 



27ri?^ ^ \5ivFkF) ^ 2^T 2' 

z, [A(r)] 

STjAjr)] ^ .„, ^^.[A(r)] , 1,^ 
_ G,n,\gn\ 1 ^ / Su^Ajr)] iz^Ajr)] 1 
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Here sum goes over the two poles. For the lowest Landau level, / = 0, 
n = Miim ( ^-^, ivF{5kF[A{r)]-kF) 1 

^ GjJ^WI /j,^ ».f(^tf[A(r)]-tr) ^ 1^ 

vr vr \ 27rT 2 

+ 27rT ^ 2V ' ^ ^ ^ 

where the shift of the Fermi momentum is given by 

(5/fcF[A(r)] = j dr./^i^^{r)^a^^^{r)^{r). (D.5) 

For T ~ Tc, we expand in A ^ T, 

n = ^^^-Im -^'( ^— ^ + -) + ^( + -) , (D.6) 

^, , Gint\qn\l^ i/dr7^V°(0^^V°(0A(r) 
A r) = Im -— — ; X 

^ + 2) + ^ + 2 V ' (^-^^ 

where ^''(x) is the derivative of the digamma function ^'{x) = ^'^^2^^ i the subleading 
term ~ A in n and the leading term ~ 1 in A vanish due to the imaginary part. We use 
that the solution of the gap equations at zero temperature is given by eq.(3.32). Therefore 
the radial profile is given by 

A(r) ~ VV'^r). (D.8) 
Substituting it into eq.(D.7), we have 

^ G,nt\qH\ 1 i (^-.,,, -ivFkF 1. ^ ,^,,,ivFkF , lA /drv/=^(V;°(r)taV0(r))2 
TT vr 2^r 2^T ^2'^^ 27rr ^ 27 R» 

(D.9) 

Simplifying the digamma functions and their derivatives, we obtain 

n=Mtanh^, (D.IO) 
2-kR 2T ' ^ ' 



^ ^ G^nt\qn\ 1 1 /c?rV=^(V^O(r)taVO(r))- 

We introduce the filling factor 



TT 2T cosh2 ^ 



(D.ll) 
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From the equation (D.IO) for the charge density, we have 

.vpkp 1 



cosh^(^) = (D-13) 



2T ' 1 



Using it in the gap equation (D.ll), we get the critical temperature for the fowest Landau 
level 

= - 4i) [ dr^g{i.\r)^a'i,\r)f. (D.14) 



For the filling factor r/-^ > 1, the critical temperature vanishes, Tc = 0, and for ry-^ < 1, 
which means either % > He or n < ric, the critical temperature grows with the magnetic 
field in the vicinity of the phase transition. The integral over the profile agrees with 
the critial temperature given in eq.(4.24). In eq.(4.24), vphi introduces the dependence 
vphi 1/ J dr g / gttip^ (r)^ il^'^ , which probably follows from a more careful definition for 
the density n in the above calculations. 
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